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Abstract 

Perspective  elongation  in  Shack-Hartmann  wavefront  sensor  (SHWFS)  spots  is  a  phe¬ 
nomenon  caused  by  an  extended  three-dimensional  scattering  beacon  in  the  meso¬ 
sphere  (80-100km).  Elongated  spots  cause  errors  in  wavefront  sensor  measurements, 
which  leads  to  poor  turbulence  compensation  and  decreased  image  resolution  of  the 
optical  system.  In  order  to  compensate  for  elongated  spots,  a  proper  beacon  model 
must  be  developed  to  simulate  the  error.  In  this  paper,  a  documented  theory  for  mod¬ 
eling  an  elongated  sodium  beacon  and  elongated  SHWFS  spots  using  sodium  layer 
“slices”  was  tested.  It  was  found  that  nine  evenly-spaced  slices  were  adequate  to  model 
the  elongated  beacon  in  the  most  stringent,  turbulence-included  case.  Furthermore 
a  bench-top  source  was  developed  and  tested  to  model  SHWFS  spot  elongation  in 
the  lab.  The  source  demonstrated  the  principle  theory,  but  requires  a  more  robust 
design  to  simulate  sodium  layer  depth.  Being  the  hrst  documented  experiment  using 
an  extended  source  on  an  adaptive  optics  (AO)  system,  it  enables  further  research  on 
the  effects  of  deep  turbulence  on  AO  systems  and  correlation  based  wavefront  sensing 
with  extended  sources. 
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SIMULATING  THE  EFFECTS  OF  AN  EXTENDED  SOURCE  ON  THE 

SHACK-HARTMANN  WAVEFRONT  SENSOR  THROUGH  TURBULENCE 

I.  Introduction 

The  Air  Force  Research  Laboratory’s  Sodium  Guidestar  Adaptive  Optics  (AO)  for 
Space  Situational  Awareness  program  (NGAS)  at  the  Starfire  Optical  Range  (SOR) 
at  Kirtland  Air  Force  Base  has  sponsored  research  on  non-uniform  spot  analysis  in 
Shack-Hartmann  wavefront  sensor  (SHWFS)  subapertures.  The  SHWFS  is  used  to 
measure  incoming  wavefront  slopes  by  measuring  the  centroid  displacement  of  focused 
spots.  When  the  spots  are  not  uniform,  errors  in  slope  measurements  occur  which 
decreases  the  resolution  of  the  optical  system.  Researching  how  to  reproduce  and 
account  for  non-uniform  spots  will  push  the  system  to  the  diffraction  limit,  and  lead 
to  higher  resolution  imagery  of  intended  targets. 

AO  is  a  technology  for  sensing  and  correcting  dynamic  wavefront  distortions  in 
real  time.  Often,  it  is  used  to  correct  for  optical  fluctuations  due  to  atmospheric 
turbulence  and  improve  imaging  system  performance.  Common  elements  in  an  AO 
system  include  a  fast-steering  mirror  (FSM),  deformable  mirror  (DM),  and  a  wave- 
front  sensor.  These  elements  are  connected  in  a  feedback,  control  loop  enabling  them 
to  communicate  and  compensate  in  real-time. 

AO  systems  require  a  light  source  or  beacon  to  measure  the  wavefront.  Ideal 
beacons  are  natural  stars  which,  when  available,  work  very  well[l].  Natural  stars 
with  the  requisite  brightness  are  not  available  in  every  part  of  the  sky,  however.  For 
the  3.5m  telescope  at  the  SOR,  a  6th  visual  magnitude  star  or  brighter  is  needed 
for  high-order  AO.  There  are  only  a  few  stars  per  degree  of  sky  coverage  with  the 
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requisite  brightness  for  AO.  At  a  wavelength  of  SOOnm,  the  probability  of  hnding  a 
bright  enough  reference  star  nearby  is  about  7  x  10“^  [2].  To  overcome  this  limitation, 
artificial  stars  created  by  lasers  mounted  on  or  near  telescope  apertures  are  used  to 
create  backscatter  in  the  atmosphere.  Typically  called  laser  guide  stars  (LGS),  these 
beacons  enable  imaging  of  object  almost  anywhere  in  sky.  Several  types  of  LGS’s 
exist,  but  this  research  focuses  on  sodium  LGS’s  since  they  are  used  at  SOR  and 
many  other  cutting-edge  astronomical  telescopes. 

In  the  past  decade,  observatories  like  SOR,  Keck,  and  others  have  begun  relying 
on  sodium  LGS’s  as  a  wavefront  reference  source[3,  4].  Sodium  LGS’s  excite  atoms 
in  the  mesosphere,  and  the  resultant  resonant  backscatter  acts  as  a  beacon  which  the 
telescope’s  AO  system  uses  to  sense  and  compensate  for  atmospheric  turbulence.  As 
an  integral  part  of  the  AO  system  at  SOR,  the  SHWFS  is  vital  in  calculating  aberrated 
wavefront  slopes.  Any  errors  occurring  within  this  device  causes  miscalculations  in 
wavefront  slopes  and,  inevitably,  lower-resolution  imagery.  This  research  focuses  on 
one  such  error:  perspective  elongation  in  SHWFS  spots.  This  research  does  not 
attempt  to  correct  this  error,  but  accurately  simulate  it  in  wave  propagation  software 
and  on  a  bench-top  laboratory  system.  Models  developed  here  be  an  important  tool 
for  researchers  who  seek  to  overcome  perspective  elongation  and  thereby  produce 
high-resolution  imagery. 

1.1  Objectives 

AO  systems  require  the  sodium  beacon  to  be  a  point  source  to  accurately  measure 
the  wavefront,  which  is  one  reason  why  natural  stars  are  desirable.  Unfortunately,  due 
to  diffraction,  the  depth  of  the  sodium  layer,  imperfect  beam  quality,  and  the  atmo¬ 
sphere,  the  beacon  scatters  from  a  large  column  and  appears  elongated  when  viewed 
by  SHWFS  subapertures.  Furthermore,  images  of  the  SHWFS  spots  are  elongated 
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differently  based  on  each  subaperture’s  viewing  angle  relative  to  the  beacon  scatter¬ 
ing  spot  [5,  6].  This  phenomenon  is  typically  referred  to  as  perspective  elongation. 
The  goal  of  this  research  is  to  explore  mitigating  solutions  for  reproducing  perspec¬ 
tive  elongation  in  SHWFS  subapertures  in  a  lab  setting.  Specihcally,  the  objectives 
include: 

1.  test  a  current,  accepted  method  used  to  simulate  a  three-dimensional  scattering 
volume  in  the  sodium  layer  with  wave  propagation  software, 

2.  quantitatively  determine  the  required  number  of  sodium  layer  slices  needed  to 
accurately  simulate  an  elongated  beacon, 

3.  create  and  implement  the  hrst  optically  designed  bench-top  method  of  produc¬ 
ing  coherent,  extended  sources, 

4.  use  the  coherent,  extended  sources  to  produce  elongation  in  a  SHWFS,  and 

5.  deliver  a  working,  scalable  model  for  simulating  perspective  elongation  in  wave 
propagation  software,  and  deliver  a  foundation  for  producing  elongation  in  an 
optics  laboratory. 

Each  of  these  goals  were  achieved  with  acceptable  results.  The  analysis  in  the  re¬ 
search  evaluated  the  number  of  sodium  layer  slices  to  sufficiently  model  an  elongated 
sodium  beacon.  The  result  is  consistent  with  Viard’s  seven-slice  model [5].  A  spin¬ 
ning  diffuser  and  adjustable  iris  were  illuminated  by  a  laser  to  create  an  incoherent, 
extended  source.  Imaging  the  source  with  the  SHWFS  produced  elongated  spots. 

1.2  Thesis  Overview 

Chapter  H  provides  an  introduction  to  Fourier  optics,  Gaussian  laser  beam  prop¬ 
agation,  sodium  layer  properties,  and  atmospheric  turbulence.  Additionally,  conven- 
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tional  AO  systems  are  discussed  including  an  in-depth  description  of  the  SHWFS 
and  the  telescope  at  SOR.  In  Chapter  III,  a  proposed  method  to  simulate  perspec¬ 
tive  elongation  using  wave  propagation  software  is  tested.  Also,  the  methodology  for 
simulating  perspective  elongation  in  a  bench-top  experiment  is  outlined.  Chapter  IV 
analyzes  the  results  of  the  software  simulations  and  lab  experiments,  and  determines 
the  most  accurate  way  to  simulate  perspective  elongation.  Finally,  Ch.  V  summarizes 
the  research  presented  here  and  highlights  key  results  and  contributions  to  science. 
Additionally,  it  discusses  ideas  for  future  work  to  further  the  research  on  perspective 
elongation. 
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II.  Background  and  Related  Research 


This  chapter  introduces  some  basic  concepts  in  the  field  of  optics.  First,  some  prin¬ 
ciples  of  Fourier  optics  are  briefly  discussed.  Next,  methods  for  modeling  atmospheric 
turbulence  are  described.  Then,  sampling  theory  for  numerical  wave  propagation  sim¬ 
ulations  is  given.  Modern  AO  systems  are  outlined,  including  a  look  at  the  system 
used  at  SOR,  and  the  theory  behind  the  Shack-Hartmann  wavefront  sensor.  Finally, 
mathematical  models  for  Gaussian  beams  are  given  along  with  characteristics  of  the 
sodium  layer  and  the  concept  of  laser  guide  stars.  This  information  provides  the 
background  necessary  to  effectively  model  perspective  elongation. 

2.1  Fourier  Optics 

This  section  introduces  some  basic  concepts  of  Fourier  optics,  which  is  a  framework 
for  describing  optical  diffraction  phenomena  using  linear  systems  theory.  Methods  of 
modeling  light  propagation  as  a  wave  are  introduced,  and  diffraction  integrals  are 
given.  Then,  the  effect  of  lenses  are  described  with  application  to  imaging  systems. 

2.1.1  Light  as  a  Wave. 

Fourier  optics  is  the  science  of  propagating  light  through  various  media.  It  differs 
from  geometrical  optics  in  that  it  models  light  as  electromagnetic  waves  instead  of 
rays,  which  accurately  accounts  for  diffraction  effects.  Figure  2  shows  the  basic  differ¬ 
ence  between  geometric  and  Fourier  optics  by  showing  how  a  plane  wave  propagates 
through  a  focusing  lens.  In  geometric  optics,  the  rays  are  focused  to  a  sharp  point  on 
the  screen,  while  in  Fourier  optics  the  wave  is  focused  to  a  Bessel  function[7]. 

When  a  light  source  is  propagated  through  space,  its  field  in  the  observation  plane 
is  found  by  numerically  evaluating  the  Fresnel  integral.  The  generalized  expression. 
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Figure  2.  In  the  geometric  optics  model  (top)  the  plane  wave  is  focused  to  a  perfect 
point  by  the  lens.  In  the  Fourier  optics  model  (bottom)  the  plane  wave  is  focused  to  a 
bessel  function  due  to  the  shape  of  the  lens  and  diffraction  effects. 


with  the  applied  paraxial  approximation,  is  given  by: 

OO  OO 

17(a;2,t/2)  =  /  j  (1) 

—  OO  — OO 


where  Xi,yi  are  coordinates  in  the  source  plane,  0:2,  t/2  are  coordinates  in  the  obser¬ 
vation  plane,  and  2;  is  the  propagation  distance [7].  By  expanding  and  rearranging  the 
terms  in  the  exponent,  Eq.  (1)  can  be  written  as: 


U{x2,y2)  = 


gjkz 

jXAz' 


J^{^2+y2) 


X 


U  (xi,  ?/i) 


— OO  —00 


(2) 


When  the  propagation  distance  Az  is  “very  far”  the  quadratic  phase  has  such  a  large 
radius  of  curvature  that  it  can  be  considered  nearly  flat  for  small  viewing  angles.  In 
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this  case,  Eq.  (2)  can  be  further  simplihed  to  the  Fraunhofer  integral  given  by: 


U  {X2,y2) 


gjkz 

jXAz 


—  OO  — OO 


(3) 


With  the  simplihcation  of  the  exponential  term  in  the  integral,  this  form  is  much 
easier  to  work  with.  To  use  the  Fraunhofer  integral,  the  propagation  distance  must 
satisfy  the  inequality: 


Az  > 


2T)2 


(4) 


where  D  is  the  source  diameter.  All  simulated  propagations  modeled  in  this  research 
require  the  Fresnel  integral  due  to  the  relatively  short  propagation  distances. 


2.1.2  Diffraction  Imaging  as  a  Linear  Systems. 

Determining  the  image  of  an  object  after  propagation  through  an  optical  system 
can  be  done  using  principles  of  linear  systems.  A  point  source  propagated  from  the 
source  plane,  through  the  optical  system,  forms  a  point  spread  function  (PSF)  at  the 
observation  plane.  The  point  source  can  be  thought  of  as  the  Greens  function  and 
the  PSF  the  spatial  impulse  response  of  the  system.  Convolving  the  PSF  with  the 
original  image  in  the  source  plane  forms  the  final  image[7].  Figure  3  shows  a  depiction 
of  this.  This  linear  system  approach  to  imaging  each  beacon  layer  is  used  in  Ch.  Ill 
to  obtain  SHWFS  lenslet  images 


2.2  Coherence  Properties  of  Light 

This  section  discusses  the  differences  between  coherent  and  incoherent  quasi- 
monochromatic  light  sources.  The  first  section  defines  Spatial  and  Temporal  co¬ 
herence  properties.  In  the  second  section,  diffusers  and  the  theory  of  transforming 
coherent  sources  into  incoherent  ones  is  outlined. 


7 


Figure  3.  Imaging  an  object  though  an  optical  system  done  by  propagating  a  point 
source  (green)  through  the  system  and  convolving  the  resultant  PSF  (red)  with  the 
original  object  (yellow)  to  get  the  final  image  (orange). 


2.2.1  Spatial  and  Temporal  Coherence. 

Coherent  properties  are  basically  measures  of  how  well  parts  of  electromagnetic 
waves  are  related  to  one  another.  Every  light  source  can  be  defined  as  spatially  coher¬ 
ent,  temporally  coherent,  both,  or  neither  (completely  incoherent).  Spatial  coherence 
is  a  measure  of  how  well  physical  areas  of  a  wavefront  are  correlated.  Temporal  co¬ 
herence  is  a  measure  of  how  quickly  two  or  more  wave  phases  drift  in  and  out  of 
sync  with  each  other.  Faster  drift  means  shorter  coherence  time,  while  slower  drift 
means  longer  coherence  time.  As  discussed  further  in  Ch.  Ill,  experiments  done  in 
this  research  only  consider  spatial  coherence  properties  of  the  source  since  temporal 
coherence  properties  don’t  have  much  effect. 


2.2.2  Diffusers. 

A  diffuser  is  an  optical  component  made  of  ground  glass  or  a  rough,  reflective 
surface.  Typically,  they  are  used  to  scatter  light  to  turn  a  spatially-coherent  source 


into  a  spatially-incoherent  one.  Diffusers  are  normally  characterized  by  grit  size  and 
pupil  size.  The  four  standard  commercial  grit  sizes,  from  coarsest  to  smoothest,  have 
designations  P120,  P220,  P600,  and  P1500  as  dehned  by  the  International  Organiza¬ 
tion  for  Standardization  (ISO).  The  grit  designations  have  average  particle  diameters 
of  125  microns,  68  microns,  25.8  microns,  and  12.6  microns,  respectively. 

When  a  quasi-monochromatic  laser  is  propagated  through  a  stationary  diffuser,  an 
unwanted  phenomenon  known  as  speckle  occurs.  Speckle  is  a  result  of  the  interference 
of  multiple  waves,  each  with  a  different,  random  phase[8].  The  addition  of  the  waves 
causes  random  constructive  and  destructive  interference  patterns  resulting  in  random 
intensities  in  the  image  plane.  Even  with  speckle,  the  spatial  coherence  of  the  laser 
is  still  intact.  A  common  way  to  remove  speckle,  and  the  spatial  coherence,  is  by 
moving  the  diffuser [9]  while  propagating  the  light  through  it.  This  basically  averages 
several  random  draws  of  the  speckle  intensity  pattern  in  the  image  plane  into  a  single, 
constant  pattern. 

2.3  Sampling  Theory 

When  performing  simulated  propagations  of  waves,  it  is  important  to  make  sure 
the  held  is  sampled  correctly.  An  example  of  unacceptable  wave  sampling  can  be 
seen  Fig.  4  where  the  ’’chirp”  function  (blue)  is  undersampled  (black).  The  “chirp” 
function  is,  coincidentally,  the  real  part  of  the  Fresnel  integral’s  quadratic  phase 
discussed  in  Sec.  2.1.  The  wave  continues  to  increase  in  frequency  but  the  sampling 
criteria  remains  constant.  In  the  beginning  the  sampling  is  hne,  but  as  time  increases 
the  wave  is  very  undersampled,  with  several  wavelengths  happening  between  each 
data  point.  This  is  called  aliasing  and  must  be  avoided  to  properly  propagate  waves. 
Typically,  one  would  apply  the  Nyquist  criterion  6  =  l/2fmax  where  fmax  is  the 
maximum  spatial  frequency  of  interest. 
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time  [s] 


Figure  4.  A  chirp  function  (blue)  with  sampling  (black).  Aliasing  is  observed  after  1 
second  because  the  frequency  of  the  signal  increases  faster  than  the  sampling. 

Following  Schmidt’s  method  for  wave  propagation  outlined  in  [10],  three  sampling 
parameters  must  be  calculated:  grid  size  N,  grid  spacing  in  the  source  plane  ()i,  and 
the  grid  spacing  in  the  observation  plane  6n-  These  values  are  dependent  on  each 
other,  and  determined  by  the  following  4  constraints: 

T  X  ^  Az-D„Si 

r.  On  ^ 

9  JV  >  ^  4-  -I- 

-  2Si  2S„  2SiS„ 

3-  (l  +  f)'5i-^<V<(l  +  f)^.  +  ^ 

where  Di  is  the  diameter  of  the  object  at  the  source  plane,  Dn  is  the  diameter  of  the 
object  at  the  observation  plane,  A  is  the  wavelength,  R  is  the  radius  of  curvature, 
z  is  the  total  propagation  distance,  and  Azi  is  the  minimum,  partial  propagation 
distance.  Since  Di,  D2,  A,  z,  and  R  are  typically  known  parameters,  the  hrst  three 
inequalities  can  be  used  to  determine  preliminary  sampling  values  for  N,  ^i,  and  Sn- 

With  the  preliminary  sampling  requirements,  the  minimum  propagation  distance 
Azi  and  the  number  of  propagations  n  can  be  determined.  While  dividing  the  total 
propagation  distance  z  into  several  smaller  propagation  distance  is  not  necessary,  it 
greatly  reduces  the  sampling  constraint  on  A^[10].  This  enables  use  of  a  relatively 
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smaller  grid  size  which  greatly  reduces  simulation  time.  Finally,  the  total  number 
of  propagations  is  simply  determined  by  n  =  ceil  {z/Azi)  +  1,  where  ceil  is  a  MAT- 
LAB  function  that  always  rounds  a  fractional  value  to  the  next  highest  integer.  It 
should  be  noted  that  this  is  the  minimum  number  of  propagations  needed,  shorter 
partial-propagation  distances  will  still  satisfy  the  fourth  inequality.  These  sampling 
constraints  are  applied  to  all  simulations  as  outlined  in  Ch.  III. 

2.4  Atmospheric  Distortion  of  Light 

The  three  atmospheric  processes  that  affect  optical  wave  propagation  are  absorp¬ 
tion,  scattering,  and  fluctuations  in  the  refractive  index[ll].  Absorption  and  scat¬ 
tering  mainly  cause  optical  wave  attenuation  while  index  of  refraction  fluctuations 
cause  irradiance  fluctuations  and  loss  of  spatial  coherence.  This  research  assumes 
absorption  and  scattering  effects  are  negligible,  but  takes  into  account  the  random, 
turbulent  refractive  index  fluctuations.  The  theory  is  described  here. 

Due  to  the  heating  and  cooling  of  the  Earth’s  surface,  large  scale  inhomogeneities 
in  air  temperature  arise  in  the  atmosphere.  These  inhomogeneities  naturally  break¬ 
down  into  smaller  and  smaller  pockets.  The  small,  randomly-evolving  pockets  are 
typically  called  eddies.  Each  eddie  in  turn  acts  as  tiny  lens  with  an  refractive  index 
that  slightly  varies  from  the  mean  value  of  1.  As  light  propagates  from  distant  objects 
in  space  through  the  atmosphere,  its  path  is  altered  and  it  accumulates  random  phase 
due  to  these  eddies.  This  leads  to  irradiance  fluctuations  and  loss  of  spatial  coherence 
in  the  wavefront.  One  can  witness  this  in  the  form  of  twinkling  stars  in  the  night 
sky.  Figure  5  shows  a  graphic  of  a  wavefront  after  it  passes  through  the  atmosphere. 
AO  systems  attempt  to  correct  these  distorted  wavefronts  to  their  ideal  form,  so  it  is 
necessary  to  effectively  model  turbulence  before  designing  an  AO  system. 
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Figure  5.  The  wavefront  from  a  distant  object  is  nearly  planar  when  it  reaches  the 
Earth’s  atmosphere.  As  the  wavefront  passes  through  the  turbulence  (gray),  it  becomes 
more  perturbed. 

2.4.1  Atmospheric  Turbulence. 

Analytically  modeling  turbulent  flow  is  challenging,  but  a  common  approach  in¬ 
volves  statistical  methods.  One  such  method  is  Kolmogorov  turbulence  theory  which 
exploits  the  fact  that  small-scale  structures  in  turbulent  air  flow  are  statistically 
isotropic  for  large  Reynolds  numbers  (Reynolds  numbers  being  the  measure  of  the 
ratio  of  inertial  forces  to  viscous  forces)  [11].  Due  to  its  sensitivity  to  temperature 
variations,  the  atmosphere’s  index  of  refraction  n(R,  t)  can  be  expressed  as 

?7,(R,  f)  =  77-0 ni(R,  t),  (5) 

where  uq  is  the  mean  value  of  the  index  of  refraction  and  ni(R,  t)  is  a  random  deviation 
from  the  mean  at  spatial  coordinates  R  and  time  t[12].  Normally,  the  mean  value  uq 
is  1  and  time  variations  are  suppressed  in  optical  wave  propagations,  reducing  Eq. 
(5)  to 

n(R)  =  1 -|- ni(R).  (6) 
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The  deviations  in  the  index  of  refraction  are  related  to  corresponding  temperature 
T(R)  and  pressure  -P(R)  fluctuations.  For  visible  and  infrared  (IR)  wavelengths  they 
are  expressed  by 


t!i(R) 


7.76  X  10“^  (1  +  7.52  x  10“V-^) 


P{R) 

r(R)’ 


For  a  wavelength  of  SOOnm,  Eq.  7  becomes: 


(7) 


n(R) 


1  +  7.99  X  10-® 


T(R)- 


(8) 


The  spatial  statistics  of  the  fluctuating  refractive  index  can  be  expressed  as  a 
spatial  power  spectral  density  (PSD).  Typically,  the  PSD  is  denoted  as  <F„(k)  where 
K,  is  the  scalar  vector  wavenumber  (rad/m)  that  describes  the  scalar  size  of  the  eddies 
in  orthogonal  components.  It  is  important  to  note  that  each  eddie  is  considered 
locally  homogenous  and  spatially  isotropic. 

Kolmogorov  theory  is  the  most  commonly  used  PSD  model  and  the  model  used 
throughout  this  study.  It  is  only  valid  for  a  limited  regime  of  turbulence  based  on 
eddie  size.  The  regime  limits  are  called  inner  scale  size  /q  and  outer  scale  size  Lq. 
Inner  scale  is  the  size  at  which  turbulence  transitions  to  laminar  flow  while  outer  scale 
is  the  size  at  which  the  eddies  are  so  large  that  they  are  no  longer  locally  homogenous 
and  isotropic.  The  expression  for  Kolmogorov  PSD  is  given  by: 


<F^(k)  =  for  ^  -C  K  < 

To  h 


(9) 


Four  other  common  PSD  models  are  the  von  Karman,  modihed  von  Karman,  Tatarskii, 
and  Hill.  These  other  methods  are  much  more  complex  and  involve  different  inner- 
and  outer-scale  factors  that  better  match  theory  with  experimental  data.  The  dif- 
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Figure  6.  The  comparison  of  different  PSD  models. 


ferent  PSD  models  are  graphically  compared  in  Fig  6.  Even  though  these  models 
are  more  accurate,  the  Kolmogorov  model  is  easier  to  work  with  and  still  provides 
sufficient  hdelity  for  this  research. 

The  parameter  in  Eq.  (9)  is  the  refractive-index  structure  parameter  measured 
in  For  horizontal  propagation  it  is  often  assumed  to  be  constant,  but  for 

vertical  propagation  it  is  a  function  altitude  h.  The  commonly  used  Hufnagel- Valley 
model  is  used  for  this  research  and  is  given  by [13]: 


Cl{h) 


0.00594  (^)hlO-ft)“exp(^) 
2.7xlO-‘»exp(j^)+Aexp(^). 


(10) 


where  v  is  the  rms  wind  speed  and  A  is  the  nominal  value  at  6*^(0) .  Setting  v  to  21 
m/s  and  A  to  1.7  x  10“^^m“^/^  results  in  the  Hufnagel- Valley  5/7  model  (H  —  V^/7), 
where  5cm  is  the  coherence  diameter  ro  and  7/irad  is  the  isoplanatic  angle  Oq.  The 
values  for  tq  and  6*0  are  calculated  at  a  wavelength  of  500nm,  although  we  adjust 
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these  numbers  for  our  wavelength,  which  is  589. 2nm.  Using  the  (if— V5/7)  model  the 
correct  tq  and  6*0  (defined  in  the  next  subsection)  are  calculated  along  with  the  log 
amplitude  variance  in  the  following  section. 

2.4.2  Atmospheric  Condition  Characterizations. 

Maxwell’s  equations  are  quantitative  expressions  that  describe  the  behavior  of 
electromagnetic  waves.  In  the  case  of  an  atmosphere-induced,  randomly-inhomogeneous, 
index  of  refraction.  Maxwell’s  equations  are  solved  by  perturbative  methods  in  con¬ 
junction  with  Green’s  functions.  One  method,  Rytov  theory,  is  typically  used  in 
conjunction  with  a  given  PSD  to  yield  statistical  moments  of  the  field  for  simple 
sources [12]  [13].  Important  moments  include  the  mean  value  of  the  field  and  the  mu¬ 
tual  coherence  function,  which  are  then  used  to  compute  properties  like  the  complex 
coherence  factor,  wave  structure  function,  phase  PSD,  and  the  mean  modulation 
transfer  function  (MTF).  The  metrics  used  in  this  research  include  the  coherence 
length  To,  the  isoplanatic  angle  6*0,  and  the  log  amplitude  variance  which  result 
from  Rytov  theory  using  the  Kolmogorov  refractive  index  PSD.  Some  of  these  param¬ 
eters  are  dependent  on  direction  of  propagation,  and  are  manipulated  for  all  relevant 
scenarios  in  this  research.  Additionally,  only  the  analytical  expressions  pertaining 
to  spherical  waves  for  these  characteristics  are  considered  since  all  propagations  are 
done  with  point  sources. 

The  coherence  length  is  a  measure  of  the  spatial  coherence  of  light.  The  larger  the 
coherence  length  the  more  coherent  the  light  is.  Figure  7  shows  a  simplihed  picture 
of  two  wavefronts  and  their  corresponding  coherence  lengths.  For  the  propagation 
of  a  point  source  from  the  ground  to  an  altitude  L  directly  overhead,  the  coherence 
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Figure  7.  The  coherence  length  estimation  of  two  aberrated  wavefronts.  The  top 
wavefront  is  less  perturbed  and,  thus,  has  a  larger  tq.  The  bottom  wavefront  is  more 
perturbed  and  has  a  smaller  tq. 


length  is  given  by[13] 

r-'/'  =  0.423fc2  J  Cl{h)  (^- j  dh,  (11) 

0 

where  k  is  the  optical  wave  nnmber.  Simply  snbstitnting  h'  =  L  —  h  into  Eq.  (11) 
gives  the  expression  for  space  to  gronnd  propagation: 


(12) 


The  typical  range  is  from  5cm  (bad  seeing)  to  20cm  (good  seeing)  for  light  at  500nm 
and  vertical  viewing  angle. 

The  isoplanatic  angle  is  a  measnre  of  how  spatially-correlated  the  atmospheric 
phase  distortion  is  for  different  object  held  angle.  A  large  angle  means  that  the 
atmosphere  is  well  correlated  for  widely-separated  object  points.  The  isoplanatic 
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Figure  8.  Two  examples  of  the  isoplanatic  angle.  In  (a)  the  entire  object  experiences  the 
same  turbulence  since  it  falls  within  6q,  but  in  (b)  the  object  experiences  uncorrelated 
turbulence  since  portions  of  it  fall  outside  Oq. 

angles  for  downward  propagation  is  given  by [13] 

L 

=  2mk^  j  Cl{h')  {L  -  h'f^  dh' .  (13) 

0 

Typical  isoplanatic  angles  range  from  5-10/trad  for  light  at  500nm  and  vertical  viewing 
angle.  Fignre  8  shows  a  graphic  of  two  different  isoplanatic  angles.  If  the  angle 
does  not  totally  encompass  the  object  as  in  Fig.  8(b)  then  different  portions  of  the 
object  lay  in  different  tnrbnlence  realizations  and  experience  nncorrelated  tnrbnlence 
effects.  This  makes  simnlating  propagations  mnch  more  difhcnlt.  Fortnnately,  all 
propagations  take  place  within  the  same  isoplanatic  angle  as  shown  in  Fig.  8(a). 

The  log  amplitnde  variance,  also  known  as  the  Rytov  nnmber  in  weak  tnrbnlence, 
measnres  the  amplitnde  flnctnations  in  propagated  light.  Typically,  this  is  very  small 
for  vertical  viewing  and,  for  point  sonrces,  is  eqnal  for  npward  and  downward  propa¬ 
gations.  The  log  amplitnde  variance  for  a  spherical  wave  is  given  by[13] 

=  0.5631A;^/®  j  Cl{h)  {L  -  hf/^dh.  (14) 

0 
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The  Rytov  number  for  weak  turbulence  usually  ranges  from  0  (no  turbulence)  to 
0.25.  While  the  Rytov  number  has  no  definitive  upper  limit,  amplitude  fluctuations 
begin  to  saturate  for  Rytov  numbers  greater  than  0.25  (strong  turbulence). 

2.5  Conventional  Adaptive  Optics 

This  section  introduces  the  basic  theory  behind  AO  systems.  It  includes  a  break¬ 
down  down  of  components  in  a  typical  AO  system  and  their  operation.  Additionally, 
it  discusses  the  system  used  on  SOR’s  3.5m  telescope. 

2.5.1  AO  Basics. 

Adaptive  optics  (AO)  is  a  technology  for  sensing  and  correcting  dynamic  wavefront 
distortions  in  real  time.  The  optical  field  of  the  wavefront  is  typically  represented  as  a 
complex  number  in  the  plane  transverse  to  propagation  with  amplitude  A  and  a  phase 
0,  and  arranged  mathematical  in  the  form  Aexp(—j(^).  As  light  propagates  through 
turbulence,  it  accumulates  an  additional  phase  (pf  Mathematically,  the  optical  held 
after  propagation  is  expressed  as 


U  =  >lexp|-j  + 


(15) 


Amplitude  huctuations  are  ignored  here  because  they  are  quite  small  for  space-to- 
ground  propagation.  Basically,  the  AO  system  senses  the  turbulent  phase  and  applies 
its  conjugate  —(pt  to  return  the  wavefront  to  its  original  form  so  that  the  compensated 
held  is  given  by: 


U  =  A  exp  [- j  (0  -1-  0t)]  X  exp  (j0t) 
=  Aexp(-j0). 


(16) 
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Figure  9.  Basic  AO  system  setup  including  the  fast-steering  mirror  (FSM),  deformable 
mirror  (DM),  and  the  wavefront  sensor  (WFS).  All  red  components  are  the  feedback, 
control  loop. 

Common  elements  in  an  AO  system  include  a  FSM,  DM,  and  a  wavefront  sensor 
(WFS).  These  elements  are  connected  in  a  feedback,  control  loop  enabling  them  to 
communicate  and  compensate  in  real-time.  Figure  9  shows  a  basic  layout  of  an  AO 
system.  An  incoming  aberrated  wavefront  is  captured  by  a  telescope  and  collimated. 
The  wavefront  sensor  detects  the  aberrations  and  commands  the  FSM  and  DM  to 
correct  for  the  distortions.  The  FSM  removes  the  overall  tilt  of  the  wavefront  and 
the  DM  corrects  for  all  higher-order  aberrations.  Tilt  typically  comprises  87%  of 
the  overall  distortion.  The  end  result  of  AO  compensation  is  typically  a  sharp,  near 
diffraction-limited  image  of  an  object.  The  scope  of  this  study  is  limited  to  wavefront 
sensor  performance,  using  a  SHWFS. 

AO  systems  can  be  operated  in  both  the  closed-loop  and  open-loop  regimes.  These 
terms  refer  to  the  configuration  of  the  control  loop.  In  closed-loop,  the  aberrated 
wavefront  reflects  off  the  DM  and  is  then  sent  through  a  beam  splitter  where  it  is 
captured  by  the  WFS  and  a  camera.  The  WFS  detects  the  wavefront  slopes  and 
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(a)  (b) 


Figure  10.  The  two  different  regimes  of  operation  for  an  AO  system,  (a)  is  the  closed- 
loops  regime  and  (b)  is  the  open-loop  regime. 

commands  the  DM  to  correct  the  wavefront.  Since  the  light  is  only  sensed  after  it 
has  been  corrected,  in  this  regime  it  is  useful  for  measuring  how  well  the  system  is 
compensating.  In  open-loop,  the  aberrated  wavefront  is  detected  before  it  has  been 
corrected  by  the  DM.  This  regime  is  useful  for  viewing  the  overall,  evolving  distortions 
of  the  wavefront.  Figure  10  shows  the  difference  between  the  two  regimes. 

2.5.2  Starfire  Optical  Range. 

SOR  is  an  observatory  located  at  Kirtland  AFB,  New  Mexico,  operated  by  the 
United  States  Air  Force  Research  Laboratory  within  its  Directed  Energy  Directorate. 
They  house  a  3.5  meter  telescope,  1.5  meter  telescope,  and  a  1.0  meter  beam  director 
with  AO  technology  designed  for  satellite  tracking.  This  requires  very  high  frame  rates 
and  a  fast  slewing  gimbal  which  is  uncommon  amongst  most  modern  observatories. 
Considered  forerunners  in  AO,  SOR  was  the  hrst  site  to  demonstrate  an  operational 
closed-loop  AO  system  utilizing  an  LGS[6][14].  Their  primary  mission  is  to  develop 
and  demonstrate  optical  wavefront  control  technologies,  and  support  field  experiments 
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Figure  11.  The  3.5  meter  telescope  at  SOR.  This  long  exposure  photo  shows  the 
exposed  telescope  aperture  and  the  sodium  LGS. 


by  others  within  the  research  community.  Figure  11  shows  a  picture  of  the  main 
telescope. 

The  AO  components  used  by  SOR  are  briefly  discussed  here.  The  3.5  meter 
Cassegrain  telescope  is  contained  within  a  retractable  dome  which  allows  the  aper¬ 
ture  to  be  complete  exposed  as  shown  in  Fig.  11.  This  enables  the  the  telescope 
to  follow  fast-slewing  objects  by  drastically  reducing  the  weight  on  the  gears.  The 
telescope’s  exit  pupil  is  imaged  onto  a  21cm  DM  containing  30  actuators  per  side.  A 
SHWFS  is  used  with  a  lenslet  array  in  the  conjugate  pupil  plane  to  sense  the  incident 
wavefront  [15].  The  beacon  for  the  AO  system  is  artificially  created  by  a  sodium  LGS 
which  can  also  be  seen  Fig.  11.  Figure  12  shows  an  image  of  the  k-Peg  binary  stars 
before  and  after  compensation.  This  image  was  captured  in  1997  and  constituted  the 
first  light  for  the  AO  system[14].  It  should  be  noted  that  the  images  are  auto-scaled, 
so  the  increased  peak  intensity  in  the  compensated  image  is  not  evident. 
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(a)  (b) 


Figure  12.  The  binary  star  system  known  as  k-Peg  as  seen  by  SOR’s  3.5  meter  telescope, 
uncompensated  image  (left)  and  compensated  image  (right).  These  images  are  auto- 
scaled,  hiding  the  increased  peak  intensity  in  the  compensation  image. 

2.6  Shack-Hartmann  Wavefront  Sensor 

This  section  explores  the  SHWFS  in  AO  systems.  First,  the  historical  develop¬ 
ment  of  the  WFS  is  discnssed.  Next,  principles  of  operation  are  explored.  Finally, 
perspective  elongation  error  is  introdnced. 

2.6.1  History  and  Development. 

The  SHWFS  is  the  and  most  common  WFS  nsed  in  AO.  The  SHWFS  is  an 
improvement  on  the  Hartmann  Screen  Test  developed  by  Johannes  Hartmann  in  the 
early  1900s.  His  simple  and  effective  methodology  involved  placing  a  hole-filled  screen 
over  the  aperture  of  the  telescope  and  studying  exposed  photographic  plates  placed  on 
either  side  of  the  focal  point  in  image  space.  The  screen  caused  light  to  pass  through 
different  entrance  pupil  locations  and  show  up  as  a  “spot  diagram  on  each  plate.  By 
comparing  the  corresponding  spot  locations  on  the  plates,  he  could  determine  where 
aberrations  resided  on  the  lens  [16]. 

After  nearly  70  years  with  little  improvement  on  the  Hartmann  Screen,  the  United 
States  Air  Force  (USAF)  became  interested  in  improving  images  of  satellites  taken 
from  Earth  in  the  1960s.  They  assigned  the  task  to  the  Optical  Science  Center  (OSC) 
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Figure  13.  The  Hartmann  screen  (left)  and  the  lenslet  array  (right).  The  lenslet  array 
uses  all  incident  photons  while  the  Hartmann  screen  only  uses  photons  incident  on  each 
hole. 


at  the  University  of  Arizona.  Dr.  Aden  Mienel  and  Dr.  Roland  Shack,  researchers 
there,  applied  the  Hartmann  Screen  Test  to  sensing  atmospheric  distortions [17]. 

Due  to  very  low  intensity  of  light  from  observed  objects  in  space.  Dr.  Shack 
replaced  the  holes  in  the  screen  with  lenses  to  focus  the  light.  He  then  optimized 
photon  use  by  placing  lenses  directly  beside  each  other  in  an  array  pattern;  a  design 
that  ensured  that  no  photon  was  wasted  by  hitting  the  space  between  the  lenses. 
Finally,  he  moved  the  lens  array  from  the  aperture  to  in  front  of  the  eyepiece.  By 
placing  a  beam  splitter  between  the  eyepiece  and  a  collimating  lens,  the  main  beam 
continued  to  the  camera  while  the  other  was  passed  to  the  array  of  lenses[17].  Figure 
13  shows  an  example  of  the  original  Hartmann  screen  and  the  modern  SHWFS  lenslet 
array. 


2.6.2  Principles  of  Operation. 

The  beauty  of  the  SHWFS  lies  in  its  simple  yet  effective  way  of  measuring  wave- 
front  slopes.  The  lenslets,  mentioned  in  the  previous  section,  are  affixed  in  front  of 
a  detector  array  as  shown  in  Fig.  14.  In  order  to  obtain  small  focal  spots  close  to 
the  diffraction  limit,  the  lenslets  are  typically  sized  so  that,  when  imaged  onto  the 
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Figure  14.  Basic  setup  of  a  Shack-Hartmann  wavefront  sensor.  Each  lenslet  focuses  a 
localized  piece  of  the  wavefront  to  a  detector  to  calculate  the  slope. 


primary  mirror,  they  are  roughly  the  size  of  the  coherence  length  tq.  Additionally, 
each  lenslet  corresponds  to  a  distinct  region  of  pixels  on  the  detector.  When  light  is 
incident  on  the  lenslet  array,  the  local  field  is  focused  onto  its  corresponding  area  of 
detector  pixels[18]. 

In  order  to  calculate  wavefront  slopes  effectively,  the  detector  array  must  be  able 
to  detect  the  centroids  of  the  focused  spots.  Two  methods  are  commonly  used:  quad 
cell  detection  and  centroid  estimation  of  light  intensity.  A  close  up  of  a  quad-cell 
sensor  from  Fig.  14  is  shown  in  Fig.  15.  Each  quadrant  is  the  size  of  a  pixel,  with 
the  sensor  centered  at  the  origin  of  the  x  —  y  plane.  The  position  of  the  focused 
beam  {T^andTy)  can  be  determined  by  comparing  the  intensity  of  light  incident  on 
the  positive  and  negative  halves  of  each  axis  as  follows: 


[C  +  D)-{A  +  B) 
""  A  +  B  +  C  +  D 
,  {A  +  D)-{B  +  C) 
^  A+B+C+D 


(17) 

(18) 
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2  pixels 

Figure  15.  An  example  of  a  quad-cell  detector. 


In  the  above  equation,  A,  B,  C,  and  D  represent  the  intensity  of  light  incident  in 
each  quadrant. 

The  other  method,  centroid  estimation,  is  performed  by  computing  the  normal¬ 
ized,  first  moment  of  the  irradiance  pattern.  Like  the  previous  method,  the  origin  is 
defined  at  the  center  of  the  subaperture  and  the  x  —  y  coordinates  of  the  focused  light 
are  calculated  by: 


rp  _  JaI  xHx,y)dxdy 

IaI 

^  yl{x,y)dxdy 

^  IaI  Hx,y)dxdy  ' 

where  A  is  the  area  of  the  subaperture  and  I{x,y)  is  the  irradiance  of  the  focused 
spot. 

Once  the  position  of  the  focused  spot  has  been  determined,  simple  geometry  is 
applied  to  determine  the  wavefront  slope  as  shown  in  Fig.  16.  In  general,  tilt  sensing 
is  a  two-dimensional  problem,  but  without  loss  of  generality,  we  can  treat  the  tilt  as 
a  one-dimensional  problem  [18].  The  slope  of  the  tilted  wavefront  with  respect  to  the 
reference  wavefront  can  be  represented  by  dW/dy.  Using  similar  triangles  as  shown 
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Figure  16.  Geometry  of  a  localized  portion  of  a  wavefront  incident  on  a  lenslet  and 
focused  on  a  detector. 


in  Fig.  16,  we  can  determine  that 


dW  _  T 
dy  R' 


(21) 


where  T  is  the  centroid  location  and  R  is  the  radins  of  cnrvatnre.  Introducing  nor¬ 
malized  pnpil  coordinates  p  =  y/r  to  Eq.  (21)  yields: 


dW  _  dW  dp 
dy  dp  dy ’ 


(22) 


where  dp/dy  is  easily  calculated  as  1/r.  Substituting  Eq.  (22)  into  Eq.  (21)  and 
rearranging  gives  an  expression  for  localized  tilt: 


dW  _  _rT 
dp  R 


(23) 


Calculating  this  for  every  lenslet  and  stitching  together  each  localized  portion  of 
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the  wavefront  effectively  reconstructs  the  original  wavefront  incident  on  the  Shack- 
Hartmann  sensor. 

2.6.3  Perspective  Elongation. 

While  several  errors  and  noise  sources  exist  within  the  SHWFS,  we  are  only  con¬ 
cerned  with  one  source  of  error:  spot  elongation.  When  an  ideal  point  source  is 
used  as  a  beacon  for  the  AO  system  in  vacuum,  the  light  is  focused  down  to  perfect, 
circular  spots  in  SHWFS  subapertures.  However,  the  spots  can  become  elongated 
due  to  diffraction,  atmospheric  effects,  imperfect  laser  quality,  and  beacon  scatter¬ 
ing  depth[2].  These  elongated  spots  lead  to  errors  in  centroid  calculations,  which  in 
turn  lead  to  incorrect  slope  measurements [6].  LGS’s  are  used  to  create  beacons  in 
the  sky  and  are  typically  launched  from  the  side  of  the  telescope.  The  amount  of 
elongation  is  proportional  to  the  separation  between  the  LGS  launching  system  and 
the  telescope  [3].  SHWFS  spots  are  elongated  differently  based  on  each  subaperture’s 
viewing  angle  relative  to  the  beacon  scattering  spot,  with  the  subapertures  furthest 
from  the  laser  showing  the  most  elongation  error [6,  5].  This  phenomenon  is  typically 
referred  to  as  perspective  elongation.  Figure  17  shows  a  simulation  of  perspective 
elongation  [4]. 

2.7  Gaussian  Beams 

Like  many  lasers,  LGS  beacons  can  be  modeled  as  a  Gaussian  beam  which  have 
several  desirable  mathematical  properties.  Gaussian  beams  remain  Gaussian  after 
propagation  through  vacuum,  although  their  beam  radius  and  wavefront  radius  of 
curvature  evolve.  Also,  beam  properties  in  the  receiver  plane  can  be  easily  calculated 
from  beam  properties  in  the  source  plane  or  vice-versa.  A  Gaussian  held  in  the  source 
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LGS  launch 
position 

Figure  17.  An  example  of  perspective  elongation.  In  this  case,  the  LGS  beacon  is 
launched  from  the  bottom  right  corner  which  causes  the  greatest  elongation  in  SHWFS 
spots  in  the  top  left  corner [4]. 
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plane  is  given  by  [12]: 


Vi  0)gaiis 


exp 


exp 


(24) 


where  Wq  is  the  source  beam  radius  and  Fq  is  the  source  radius  of  curvature. 

For  any  propagation  distance  normal  to  the  x  —  y  plane,  the  Gaussian  held  and 
intensity  are  given  by  [12]: 


U{x,y,z) 

I{x,y,z) 


gaus 


gaus 


1  /  \ 
00  +  j Ao 

1  /  2r2\ 


(25) 

(26) 


where  ©o  is  the  amplitude  change  due  to  focusing,  Aq  is  the  amplitude  change  due 
to  diffraction,  W  is  the  beam  radius  at  z,  and  F  is  the  radius  of  curvature  at  z.  As 
stated  above,  these  beam  parameters  can  be  calculated  from  source  parameters  as 
follows: 


©0  {z) 
Ao  (z) 
W{z) 
F{z) 


2z 

W„^Ql  +  hl 
fo  (eg  +  Ag)  (90  -  1) 
®0  +  Aq  ~  ©0 


(27) 

(28) 

(29) 

(30) 


While  the  above  theory  applies  to  Gaussian  beam  propagation  in  vacuum,  an 
expression  for  the  mean  turbulence-degraded  irradiance  of  an  initially  Gaussian  beam 
has  been  developed.  Glosed-form  expressions  for  Gaussian  beams  through  turbulence 
are  involved  to  derive,  but  one  common  expression  is  mean  intensity  at  any  given 
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propagation  distance  given  by [12]: 


(/(r,.))  =  ^exp 


2r2 

'm 


exp  [2cr^(r,2;)  -  T]  , 


(31) 


where  describes  the  atmospherically  indnced  change  the  transverse  direction,  and 
T  describes  the  change  in  the  on-axis  mean  irradiance  at  the  receiver  plane  cansed  by 
turbulence.  Both  terms  are  proportional  to  the  reverse-path  tq  value  from  Eq.(12)  as 
evident  in  the  expressions: 


T 


3.62fc^/6^5/6^5/6 


(32) 

(33) 


With  a  minor  rearrangement  of  Eq.  (12),  Eqs.  (32)  and  (33)  can  be  rewritten  in 
closed  form  as: 


T 


8.30 

9.98 


AL\ 

krl) 

\  5/6 


5/6 


(34) 

(35) 


It  should  be  noted  that  Eq.  (31)  is  the  ensemble  average  irradiance  (also  called 
long-exposure)  that  assumes  the  beam  is  not  tracked. 


2.8  Laser  Guide  Stars  and  Sodium  Layer  Properties 

As  shown  in  Fig.  9,  the  AO  system  requires  a  light  source  or  beacon  to  gather 
wavefront  data.  And,  as  mentioned  in  the  introduction,  observatories  like  SOR,  Keck, 
and  others  have  begun  relying  on  sodium  LGS’s  as  a  wavefront  reference  source[3,  4]. 
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Figure  18.  An  elongated  beacon  caused  by  the  depth  h^a  of  the  sodium  layer  results 
in  elongated  spots  in  SHWFS  subapertures. 

Figure  11  shows  the  beam  of  the  sodium  LGS  launched  at  the  SOR.  At  a  wavelength 
of  589. 2nm  and  separated  a  distance  a  from  the  center  of  a  telescope  aperture  H, 
sodium  LGS’s  excite  atoms  in  the  mesosphere  as  shown  in  Fig.  18. 

Using  the  sodium  layer  to  create  artificial  beacons  has  several  advantages  and 
disadvantages.  To  fully  understand  them,  sodium  layer  properties  are  described  here. 
The  mesosphere  contains  a  buildup  of  neutral  sodium  atoms  located  at  a  mean  height 
FTjva  of  90km.  The  depth  of  the  sodium  layer  can  vary  from  5-20km  depending 
on  location,  season,  and  even  time  of  day [19,  20].  While  the  origins  of  sodium  atom 
buildup  in  the  mesosphere  is  unknown,  it  has  been  hypothesized  that  it  was  formed 
from  the  ablation  of  meteors[21].  Studies  have  shown  that  the  distribution  of  sodium 
atoms  is  Gaussian  as  shown  in  Fig.  19[6,  21].  At  different  distances  Hi  from  the  the 
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100km  90km  80km 


Figure  19.  Model  column  density  of  sodium  atoms  in  the  mesosphere  at  a  mean  height 
H^a  of  90km. 


mean,  the  column  density  can  be  approximated  by 


/  {Hl)  =  exp 


(Hm  - 

m 


(36) 


where  Wi  from  the  full- width,  half-max  as  5km/\/ln  2. 

When  the  sodium  atoms  are  excited  by  the  laser,  they  emit  light  becoming  a 
beacon  for  the  AO  system.  These  beacons  can  be  placed  anywhere  in  the  sodium 
layer,  theoretically  enabling  the  system  to  image  anything  in  the  sky,  including  fast- 
moving  objects.  Additionally,  the  beacons  are  brighter  than  natural  stars  which  gives 
the  AO  system  an  ample  amount  of  photons  to  work  with.  The  problem  comes  with 
the  depth  of  the  sodium  layer.  As  mentioned  before,  an  ideal  beacon  is  a  point  source, 
but  the  depth  forms  an  oblong  three-dimensional  scattering  volume  when  excited  by 
the  laser  leading  to  one  of  the  causes  of  perspective  elongation  discussed  in  Sec.  2.6. 
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2.9  Imagery  Comparison  Techniques 


Perspective  elongation  is  simulated  multiple  times  using  different  methods  of  mod¬ 
eling  an  elongated  beacon  in  the  sodium  layer.  To  determine  the  optimal  elongated 
beacon  model,  it  was  necessary  to  compare  the  elongated  spot  images  in  SHWFS  sub¬ 
apertures.  The  two  methods  chosen  in  this  research  were  root  mean  squared  (RMS) 
error  and  correlation.  The  following  subsections  explain  both  methods  in  more  detail. 

2.9.1  Root  Mean  Square  Error. 

The  RMS  error  is  typically  used  to  measure  the  differences  between  simulated 
values  from  a  model  and  actual  values  from  the  thing  being  modeled.  This  method 
determines  the  model’s  precision.  The  mathematical  definition  is  given  by 


(37) 


RMSiiff 


where  7^  is  an  arbitrary  symbol  for  a  data  set  or  value  and  (...)  is  the  ensemble 
average.  A  slight  variation  on  the  RMS  error  is  the  RMS  percent  difference.  This 
method  is  essential  the  same  as  the  RMS  error  calculation  except  it  gives  the  percent 
difference.  The  mathematical  definition  is 


(38) 


2.9.2  Correlation. 

Another  useful  comparison  technique  is  called  correlation.  Correlation  is  a  statistic 
measuring  the  dependence  between  two  or  more  data  sets.  One  way  to  measure 
correlation  is  by  calculating  the  correlation  coefficient  p.  This  is  a  number  that  ranges 
from  from  -1  to  1  where  -1  is  completely  anticorrelated,  0  is  completely  uncorrelated. 
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and  1  is  completely  correlated.  Values  close  to  -1  or  1  are  considered  highly  correlated 
and  thus  very  dependant,  while  values  close  to  0  are  considered  highly  uncorrelated 
and  thus  independent.  The  mathematical  expression  for  the  correlation  coefficient  is 
given  by: 

cov{A,B)  _  {AB)  -  {A){B) 
ctactb  o-actb 

where  A  and  B  are  data  sets,  (...)  is  the  ensemble  average,  and  a  is  the  standard 

deviation.  Since  the  dehnition  manipulates  different  statistical  moments,  it  is  very 

useful  for  comparing  images  formed  by  the  data  instead  of  the  raw  values  of  the  data 

themselves. 

2.10  Chapter  Summary 

In  summary,  AO  systems  use  beacons  to  collect  data  on  distorted  wavefronts. 
The  distortions  are  caused  by  inherent  turbulence  in  Earth’s  atmosphere.  Although 
ideal  beacons,  natural  stars  of  the  required  brightness  are  rare,  so  Gaussian  lasers 
are  used  to  create  artihcial  beacons  in  the  sodium  layer  of  the  atmosphere.  This 
enables  optical  systems,  like  the  observatory  at  SOR,  to  image  almost  anything  in 
the  sky.  Unfortunately,  the  artificial  beacons  are  distorted  due  to  the  depth  of  the 
sodium  layer  and  the  separation  from  the  telescope  aperture.  A  distorted  beacon 
causes  non-uniform  spots  in  SHWFS  subapertures.  This  leads  to  an  inaccuracy  of 
wavefront  slope  measurements  and  a  lost  in  the  image  resolution  of  the  optical  system. 
Chapter  III  discusses  methods  used  to  accurately  simulate  the  distorted  beacon  and 
form  elongated  spots  in  SHWFS  subapertures. 


Pa,b  =  corr(A,  B) 
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III.  Research  Methodology 


This  chapter  explains  describes  the  methodology  for  both  simulations  and  exper¬ 
iments  used  in  this  research.  The  goal  of  this  research  was  to  develop  models  that 
accurately  simulated  perspective  elongation  in  SHWFS  subapertures,  and  was  done 
in  two  parts.  The  first  part  used  MATLAB  to  model  the  elongation  by  testing  a 
method  described  in  current  literature.  The  simulation  was  used  to  determine  how 
accurate  the  documented  method  was.  The  second  part  used  a  special  source  in  an 
optics  laboratory  to  recreate  the  elongation  in  a  real  SHWFS.  The  experimental  result 
was  then  compared  to  the  MATLAB  simulation  for  validation. 

3.1  Computer  Simulation  Environment 

To  study  perspective  elongation  in  SHWFS  subapertures,  it  was  necessary  to  sim¬ 
ulate  an  elongated  beacon  in  the  sodium  layer.  This  was  done  by  simulating  the 
propagation  of  a  Gaussian  beam  through  the  atmosphere,  from  ground  to  sky.  Con¬ 
trolling  how  the  beam  interacted  with  the  sodium  layer  created  the  elongated  beacon. 
Once  the  beacon  was  manipulated  to  resemble  the  beacon  shape  observed  at  SOR, 
it  was  propagated  back  through  the  atmosphere  and  the  SHWFS.  For  simplicity,  a 
telescope  was  omitted  in  all  simulations  and  the  SHWFS  lenslet  array  was  placed  at 
the  primary  aperture  plane  instead.  The  telescope  was  not  necessary  since  it  only 
demagnifies  the  image  and  has  no  direct  effect  on  perspective  elongation.  Addition¬ 
ally,  all  simulations  were  done  in  open-loop  as  described  in  Sec.  2.5  to  show  the 
evolution  of  the  elongated  spots.  Figure  20  shows  a  basic  layout  of  the  simulation 
geometry  including  the  parameters.  It  should  be  noted  that  a  2m  aperture  separated 
4m  from  the  LGS  launching  system  was  used  instead  of  the  SOR  specifications  of 
a  3.5m  aperture  separated  2.5m  from  the  LGS  launching  system.  Since  the  total 
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Figure  20.  The  basic  setup  of  the  software  model  geometry.  The  telescope  is  replaced 
with  SHWFS  lenslets  to  simplify  the  setup. 

distance  for  both  scenarios  is  the  same  (6m),  our  setup  basically  recreates  the  same 
geometry  used  at  SOR.  The  smaller  aperture  greatly  reduced  simulation  times  while 
keeping  the  subaperture  geometry  the  same. 

Elongated  beacons  in  the  sodium  layer  are  difficult  to  simulate  because  they  are 
three-dimensional  scattering  volumes.  Research  shows  that  they  can  be  simulated 
by  dividing  the  sodium  layer  into  multiple,  evenly-spaced  slices  and  propagating  a 
Gaussian  beam  from  a  fixed  location  on  the  ground  to  each  slice  [5].  This  creates 
a  different  intensity  pattern  in  each  individual  slice.  In  past  studies  seven  slices  in 
the  sodium  layer  were  used  without  verification  to  create  these  scattering  volumes[5]. 
Figure  21  shows  the  division  of  the  sodium  layer  into  seven  slices. 

While  using  seven  slices  does  form  the  elongated  beacon  and  recreates  the  per¬ 
spective  elongation  in  SHWFS  subapertures,  the  accuracy  of  this  configuration  has 
not  been  reported.  To  test  this  method,  the  elongated  beacon  was  modeled  19  times. 
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Beacon  Spot  Area 


Figure  21.  The  sodium  layer  divided  into  discrete  slices  between  80km  and  100km. 

each  time  using  an  increasing,  odd  number  of  evenly-spaced  sodium  layer  slices  from 
3  to  39.  Each  elongated  beacon  was  then  imaged  by  a  SHWFS,  resulting  in  elongated 
subaperture  images.  Comparing  the  elongated  spots  revealed  the  minimum  number 
of  sodium  layer  slices  needed. 

To  insure  the  perspective  elongation  was  recreated  properly,  several  pieces  were 
individually  simulated  and  verified.  Since  the  WaveProp  software  was  used  for  nearly 
all  facets  of  software  simulation,  typical  verification  entailed  performing  the  desired 
propagation/modeling  in  WaveProp  and  comparing  the  results  to  theory.  The  follow¬ 
ing  subsections  provide  the  verification  of  the  following  parameters:  grid  sampling, 
atmospheric  turbulence  characteristics,  and  Gaussian  beam  propagation. 
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Table  1.  Sampling  parameters  for  all  wave  propagations. 


Parameters 

Upward 

Downward 

Grid  size  N 

1024 

1024 

Grid  spacing  (source)  (5i  [mm] 

4 

4 

Grid  spacing  (observation)  5n  [mm] 

10 

10 

Max  propagation  distance  Azmax  [km] 

27.8 

27.8 

Number  of  propagations  n 

10 

10 

3.1.1  Grid  Sampling. 

To  properly  simulate  the  propagation  of  light  on  a  hnite  grid,  the  sampling  pa¬ 
rameters  were  carefully  determined  from  the  method  outlined  in  Sec.  2.3.  For  the 
propagation  of  the  Gaussian  beam  from  ground-to-sky,  the  source  diameter  Di  was 
set  to  8cm,  the  diameter  of  the  object  in  the  observation  plane  Dn  was  set  to  1.25m, 
and  the  propagation  distance  z  was  set  to  90km.  These  lengths  were  similar  to  ac¬ 
tual  specihcations  of  the  LGS  AO  system  used  at  SOR.  To  satisfy  all  four  sampling 
inequalities,  they  were  plotted  on  the  same  chart  and  analyzed.  As  shown  in  Fig. 
22  the  contour  lines  represent  different  grid  sizes  N  and  the  dashed  line  represents 
constraint  1.  Gonstraints  2  and  3  are  not  shown  because  they  plotted  well  out  of  the 
range  and  essentially  do  not  impact  the  sampling  analysis.  The  sampling  parameters 
from  ground-to-sky  are  more  stringent  than  those  from  sky-to-ground,  so  the  same 
parameters  were  used  for  all  propagations.  Table  1  lists  the  parameters  used  in  the 
simulations. 

3.1.2  Atmospheric  Characteristics  and  Verification. 

To  accurately  simulate  perspective  elongation,  an  atmospheric  model  was  devel¬ 
oped.  As  outlined  in  Sec.  2.4,  the  atmosphere  can  be  dehned  with  quantitative 
parameters  such  as  the  coherence  length  ro,  isoplanatic  angle  and  the  Rytov  num¬ 
ber  (T^.  Using  the  equations  in  Sec.  2.4  the  atmospheric  parameters  for  this  study  are 
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5j  [m] 


Figure  22.  A  plot  of  the  sampling  inequalities.  The  contour  lines  are  different  grid 
sizes  N  and  the  dashed  line  is  the  maximum  boundary  for  the  grid  spacing  parameters 
Si  and  Sn-  The  data  point  shows  the  chosen  parameters. 


39 


Table  2.  Theoretical  atmospheric  parameters  for  vertical  upward  and  downward  prop¬ 
agation  paths. 


Metric 

upward 

downward 

ro  [m] 

2.40 

0.0619 

6*0  [/irad] 

— - 

8.402 

0.045 

0.045 

shown  in  Table  2  for  both  the  gronnd-to-sky  (npward)  and  sky-to-gronnd  (downward) 
propagations.  All  characteristics  were  calculated  were  point  sources.  A  wavelength 
of  589.2nni  and  a  propagation  distance  L  =  QOkm  was  used  in  all  calculations.  As 
evident  in  the  table,  the  parameters  fall  within  their  expected  ranges.  Since  all  val¬ 
ues  were  calculated  for  vertical  propagation  paths,  the  overall  turbulence  is  relatively 
weak,  with  a  Rytov  number  of  only  0.045. 

The  atmospheric  parameters  were  then  used  to  program  a  three-dimensional  tur¬ 
bulence  model  in  WaveProp.  A  three-dimensional  model  more  accurately  simulates 
the  effects  of  turbulence  and,  in  this  case,  consisted  of  10  phase  screens  placed  along 
the  propagation  path.  To  verify  that  the  model  accurately  simulated  Kolmogorov 
turbulence,  the  coherence  length  was  tested  for  a  single  phase  screen.  Finally,  tests 
to  verify  the  coherence  length  and  Rytov  number  for  the  entire  three-dimensional 
model  were  performed. 

To  verify  tq  for  a  single  phase  screen,  the  phase  structure  function  was  com¬ 
puted  and  compared  to  theory.  For  the  Kolmogorov  spectrum,  the  phase  structure 
function  is  given  by  [10]: 

D^(Ar)  =6.88  f  —  j  ,  (40) 

where  Ar  is  the  spatial  separation  between  points.  Using  the  coherence  length  tq  for 
downward  propagation  in  Table  2,  phase  screens,  such  as  the  one  shown  shown  in 
Fig.  23,  were  generated  in  WaveProp.  Because  of  the  random  nature  of  phase  screen 
generation,  100  phase  screens  and  their  corresponding  phase  structure  functions  were 
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Figure  23.  A  randomly-generated  Kolmogorov  phase  screen. 

computed  and  averaged.  Figure  24  shows  a  plot  of  the  theoretical  and  simulated 
D^{Ar).  As  evident  in  the  hgure,  the  functions  match  up  as  well  as  the  simulation 
method  allows. 

Similarly  to  the  single  phase  screen,  the  3D  turbulence  model  was  verihed  with 
the  same  phase  structure  function  ^^(Ar)  in  Eq.  (40).  Using  the  same  ro  value,  100 
phase  structure  functions  were  simulated  and  averaged  in  MATLAB.  Figure  25  shows 
a  plot  of  the  theoretical  and  simulated  D^{Ar).  The  functions  match  up  well,  thus 
verifying  that  the  turbulence  model  is  behaving  as  anticipated. 

To  verify  the  cr^,  the  log-amplitude  PDF  was  calculated  first.  An  on-axis  point 
source  located  at  a  height  of  90km  from  the  aperture  plane  ws  propagated  through  a 
single  Kolmogrov  phase  screen  and  vacuum.  The  on-axis  log  amplitude  y  was  then 
calculated  by [8,  12]: 

^  (^) ' 
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Figure  24.  The  theoretical  (blue)  and  the  simulated  (black)  structure  functions  for 
a  single  Kolmogorov  phase  screen.  The  simulated  structure  function  was  found  by 
averaging  100  randomly  generated  structure  functions. 


Figure  25.  The  theoretical  (blue)  and  the  simulated  (black)  structure  functions  for  the 
three-dimensional  (10  screen)  Kolmogorov  turbulence  model.  The  simulated  structure 
function  was  found  by  averaging  100  randomly  generated  structure  functions. 
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Figure  26.  Comparison  of  the  on-axis  log  amplitude  with  a  Gaussian  PDF. 


where  \Ut\  is  the  amplitude  of  the  held  through  turbuleuce  aud  |f/o|  is  the  amplitude 
of  the  held  through  vacuum.  Due  to  the  raudom  uature  of  the  turbuleuce  model,  1000 
realizatious  of  the  propagatiou  were  executed,  aud  the  log  amplitude  was  calculated 
for  each.  Ouly  the  ou-axis  grid  poiuts  were  used  iu  each  calculatiou.  Results  were 
plotted  iu  a  histogram  aud  compared  to  a  Gaussiau  PDF  as  showu  iu  Fig.  26.  To 
further  verify  the  Rytov  uumber,  the  RMS  error  was  tested  for  several  Rytov  uumbers 
betweeu  0  aud  0.1  usiug  the  results  from  the  histogram.  Figure  27  shows  a  plot  of 
the  RMS  error.  The  resultaut  Rytov  uumber  is  close  to  the  theoretical  value  iu  Table 
2.  With  a  valid  turbuleuce  model,  propagatious  of  light  through  the  atmosphere 
will  behave  as  expected  aud  cau  be  verified  with  theory  as  discussed  iu  the  uext 
subsectiou. 
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Figure  27.  The  mean  squared  error  plot  to  locate  the  simulated  Rytov  value. 


3.1.3  LGS  Beam  System  Parameters. 

The  laser  used  in  the  LGS  system  was  modeled  as  a  Gaussian  beam.  Using  the 
specihcation  from  SOR,  the  source  diameter  Wq  was  set  to  8cm  and  the  source  radius 
of  curvature  Fq  was  set  to  90km  [19,  20].  These  variables  were  kept  constant  when 
calculating  the  source  amplitude  parameters  due  to  focus  ©o  and  diffraction  Aq  for 
different  propagation  distances  z.  The  beam  parameters  W  and  F  in  the  observation 
plane  were  then  calculated  from  the  source  parameters.  Propagation  distances  were 
set  by  the  number  of  slices  used  to  divide  the  sodium  layer.  Table  3  shows  data  for 
all  parameters  used  in  the  7  slice  scenario.  It  should  be  noted  that  the  Gaussian 
laser  was  not  focused  to  the  mean  of  the  sodium  layer  (90km).  Figure  28  shows  the 
spot  size  trend  for  vacuum  propagation  of  the  beam  from  ground  to  sky.  The  graph 
shows  that  the  beam  focuses  around  10km,  well  before  the  mean  of  the  sodium  layer. 
However,  the  beam  waist  only  increases  approximately  5cm  over  the  entire  depth 
of  the  sodium  layer  (20km).  The  code  used  to  model  Gaussian  beam  parameters  is 
shown  in  Listing  A.l  of  Appendix  A. 

The  source  beam  parameters  were  then  programmed  into  WaveProp  and  prop¬ 
agated  to  various  slices  (distance)  in  the  sodium  layer.  This  effectively  formed  the 
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Table  3.  Theoretical  data  for  Gaussian  beam  vacuum  propagation  from  ground  to  sky. 
This  is  for  the  scenario  of  dividing  the  sodium  layer  into  7  slices. 


propagation  distance  2;  [km] 


80 

83.333 

86.667 

90 

93.333 

96.667 

100 

Wq  [cm] 

8 

8 

8 

8 

8 

8 

8 

Fo  [km] 

90 

90 

90 

90 

90 

90 

90 

Aq 

2.4909 

2.5397 

2.5886 

2.6374 

2.6862 

2.7351 

2.7839 

00 

0.0556 

0.037 

0.0185 

0 

-0.0185 

-0.037 

-0.0556 

W  [cm] 

19.93 

20.32 

20.71 

21.10 

21.49 

21.88 

22.28 

F  [km] 

-85.768 

-87.176 

-88.578 

-90 

-91.432 

-92.874 

-94.324 

propagation  distance  [km] 


Figure  28.  A  graph  showing  the  trend  of  the  spot  size  from  ground  to  the  outer  edge  of 
the  sodium  layer  (100km).  The  waist  forms  around  10km,  but  the  beam  size  increases 
very  slowly  through  the  entire  propagation. 
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Figure  29.  The  center  cutaway  of  the  Gaussian  beam  vacuum  propagation  to  7  slices 
of  the  sodium  layer.  The  simulation  result  (blue)  match  perfectly  with  the  theoretical 
expressions  (black). 


elongated  beacon  model.  To  verify  the  beams  propagated  correctly,  the  t/  =  0  cut¬ 
away  of  the  intensity  pattern  at  each  slice  was  compared  to  theory.  This  was  done 
in  vacuum  propagation  for  the  7  slice  scenario  and  shown  in  Fig  29.  As  evident  in 
the  hgure,  the  models  match  perfectly  with  the  theoretical  expressions.  Listing  A. 2 
in  Appendix  A  shows  the  code  used  to  propagate  and  verify  Gaussian  beams. 

For  the  turbulence  scenario,  the  beams  were  programmed  into  WaveProp  and 
propagated  through  the  turbulence  model.  Due  to  the  randomness  of  the  turbulence 
model,  each  propagation  yielded  different  results.  To  verify  that  the  propagations 
were  working  properly,  10000  propagations  were  performed  and  averaged.  This  was 
then  compared  to  the  turbulent  degraded  irradiance  pattern  outlined  in  Sec.  2.7.  For 
simplicity,  the  propagation  distance  was  set  to  90km  for  each  iteration.  Figure  30 
shows  the  comparison  of  the  WaveProp  model  and  theory.  The  models  match  as  well 
as  the  simulation  allowed.  If  the  number  of  propagations  were  increased  to  inhnity, 
the  plots  would  match  up  perfectly.  Listing  A. 3  in  Appendix  A  shows  the  code  used 
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Figure  30.  The  center  cutaway  of  the  Gaussian  beam  turbulence  propagation  to  the 
mean  of  the  sodium  layer  (90km).  The  simulation  result  (blue)  has  the  same  trend  as 
the  theoretical  expression  (black). 

to  propagate  and  verify  turbulence-degraded  Gaussian  beams.  Now  that  our  method 
to  create  an  elongated  beacon  model  has  been  verihed,  the  next  step  is  to  view  the 
beacon  with  the  SHWFS  lenslets  to  show  perspective  elongation.  The  next  subsection 
explains  how  this  was  accomplished. 


3.1.4  Simulating  Perspective  Elongation. 

To  simulate  perspective  elongation  in  SHWFS  subapertures,  the  hrst  step  was  to 
propagate  the  elongated  beacon  back  through  the  atmosphere  (sky  to  ground)  to  the 
telescope  aperture  plane.  As  a  reminder,  an  actual  telescope  was  not  used  since  it 
has  no  direct  effect  on  the  elongation  phenomenon.  Instead,  the  SHWFS  lenlets  were 
placed  in  the  primary  telescope  aperture  plane. 

The  Fourier  optics  principles  in  Sec.  2.1  were  exploited  for  all  downward  propaga¬ 
tions.  Instead  of  a  complex  shape,  point  sources  were  propagated  through  the  entire 
optical  system  to  obtain  the  PSF’s  in  each  subaperture.  They  were  then  convolved 
and  averaged  with  the  original  objects  to  obtain  the  hnal  image.  The  “objects”  in 
this  case  were  the  irradiance  patterns  formed  by  the  Gaussian  beam  propagation  to 


47 


y  [m]  y  [m] 


Figure  31.  Gaussian  beam  irradiance  pattern  at  central  slice  of  the  sodium  layer  (90km). 
Propagated  through  vacuum  (left)  and  through  turbulence  (right). 

different  slices  (distances)  in  the  sodium  layer  (as  discussed  in  the  previous  subsec¬ 
tion).  As  an  example,  Fig.  31  shows  the  irradiance  patterns  through  vacuum  and 
turbulence  of  the  center  slice  of  the  sodium  layer  (90km).  It  should  be  noted  that  the 
same  turbulence  draw  was  used  in  all  propagations.  Also,  the  tilt  was  removed  from 
the  turbulence  beacon  in  the  downward  propagations.  This  can  be  seen  in  Fig.  31, 
where  peak  of  the  perturbed  irradiance  pattern  is  centered  at  the  origin. 

The  number  of  point  sources  used  corresponded  to  the  number  of  slices  the  sodium 
layer  was  divided  into.  For  example,  in  the  seven  slice  scenario,  a  point  source  was 
propagated  from  each  slice  through  the  turbulent  path  and  the  optical  system.  The 
separation  distance  a  between  the  LGS  launching  system  and  the  SHWFS  lenslets 
was  modeled  by  shifting  the  point  sources  in  each  slice  off  the  optical  axis.  Fixing 
the  mean  slice  at  90km  aligned  with  the  optical  axis,  point  sources  were  shifted 
in  the  x  and  y  directions  by: 
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Figure  32.  The  point  source  coordinates  determined  by  basic  trigonometry.  This  figure 
only  shows  the  coordinates  for  the  point  source  in  the  top,  middle,  and  bottom  slices 
for  simplicity. 


Table  4.  An  example  of  the  point  source  coordinates  for  the  seven-slice  scenario. 


Slice 


X 


m 


y 


m 


2;  [km] 


1 

-0.444 

-0.444 

80 

2 

-0.296 

-0.296 

83.333 

3 

-0.148 

-0.148 

86.667 

4 

0 

0 

90 

5 

0.148 

0.148 

93.333 

6 

0.296 

0.296 

96.667 

7 

0.444 

0.444 

100 

where  Hs  is  the  height  of  a  slice.  These  equations  were  derived  from  the  basic 
geometric  setup  of  the  simulation.  Figure  32  shows  the  geometry  for  point  source 
separation,  and  Table  4  lists  the  point  source  coordinates  for  the  seven-slice  scenario. 


After  propagating  the  point  sources  through  the  optical  system,  PSF’s  were  formed 
in  each  SHWFS  subaperture.  Before  the  PSF’s  were  convolved  with  the  original  ob¬ 
ject,  the  object’s  phase  was  collimated  by  removing  the  spherical  vacuum  phase  of 
the  center  slice  from  every  other  slice.  This  was  done  to  make  the  SHWFS  lenslets 
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Figure  33.  Elongated  SHWFS  spots  formed  from  the  seven  slice  model,  (top)  are  spots 
through  vacuum  and  (bottom)  are  spots  through  turbulence.  This  is  a  small  sample 
(16  subapertures)  of  the  entire  SHWFS  array  (400  subapertures). 

focused  to  the  mean  distance  of  the  sodium  layer  (90km)  since  a  telescope  was  not 
used.  Finally,  after  the  convolution  was  performed  and  averaged  for  each  PSF  and 
object  from  different  sodium  layer  slices,  the  perspective  elongation  emerged.  Propa¬ 
gations  were  done  through  both  vacuum  and  turbulence  conditions,  and  an  example 
of  the  elongated  spots  in  the  SHWFS  is  displayed  in  Fig.  33.  This  figure  shows  the 
spots  formed  by  the  seven  slice  model,  and  only  a  small  portion  of  the  entire  SHWFS 
array  is  displayed  (16  subapertures). 

The  entire  process  was  repeated  to  simulate  perspective  elongation  using  different 
numbers  of  sodium  layer  slices.  The  elongated  spot  images  were  compared  by  using 
RMS  percent  error  and  correlation  coefficient.  Listing  A. 4  in  Appendix  A  displays 
the  code  used  to  simulate  perspective  elongation.  The  results  are  analyzed  in  Ch. 
IV.  With  a  valid  model  for  simulating  perspective  elongation  in  WaveProp,  the  next 
step  was  to  experimentally  create  elongated  spots  in  a  SHWFS.  The  experimental 
methodology  is  discussed  in  the  next  section. 
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3.2  Laboratory  Description 


The  purpose  of  this  section  is  to  introduce  the  AO  lab  bench-top  setup  used 
throughout  this  study.  First,  the  operation  of  the  AO  system  is  discussed.  Then, 
a  ZEMAX  model  of  the  bench-top  was  then  developed  to  determine  the  theoretical 
magnification  and  other  characteristics  of  the  system. 

3.2.1  AFIT  AO  Lab. 

The  AO  lab  at  AFIT  contains  all  the  components  necessary  to  operate  an  AO 
system  in  the  open-  and  closed-loop  regimes.  The  system  is  capable  of  performing 
in  moderate  to  strong  turbulence,  with  Rytov  numbers  around  1.  A  diagram  of 
the  optical  setup  is  shown  in  Fig.  34,  and  the  components  are  listed  in  Table  5. 
The  phase  wheels  and  lenses  L1-L7  make  up  the  atmospheric  turbulence  simulator 
(ATS).  The  ATS  is  capable  of  producing  a  broad  range  of  dynamic  and  static,  low 
altitude  and  high  altitude  turbulence  scenarios.  A  fiber-coupled  laser  diode,  creating 
a  near-point,  is  propagated  through  the  ATS  and  reflected  off  the  FSM.  The  FSM 
applies  tip/tilt  corrections  to  the  collimated  output  beam  based  on  commands  made 
using  centroid  measurements  from  the  tracker  camera.  After  the  FSM,  the  collimated 
beam  enters  a  4-f  lens  system  composed  of  lenses  L8  and  L9.  The  4-f  system  creates 
a  conjugate  pupil  plane  between  the  exit  pupil  of  the  ATS  and  the  surface  of  the  DM. 
The  ratio  of  the  focal  lengths  of  lenses  L8  and  L9  applies  proper  magnification  to 
shrink  the  beam  to  match  the  DM  size.  Reflection  off  the  DM  then  applies  higher- 
order  aberration  corrections  based  on  commands  made  from  gradient  measurements 
using  the  SHWFS.  The  collimated  light  is  then  propagated  through  another  4-f  system 
(with  unit  magnification)  and  a  series  of  beam  splitters  and  folding  mirrors  until  it 
is  captured  by  the  tracker  camera  and  the  SHWFS.  All  optical  components  past  the 
ATS  comprise  the  AO  system. 
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Figure  34.  Physical  layout  of  the  turbulence  simulator  and  AO  system. 


Table  5.  Lab  bench  components 


Symbol 

component 

LD 

laser  diode 

Ll-Lll 

lenses 

PHI  1  and  PHI  2 

phase  wheels 

Ml,  M3,  and  M4 

fold  mirrors 

BSl  and  BS2 

beam  splitters 

FSM 

fast  steering  mirror 

DM 

deformable  mirror 

Tracker 

tracking  camera 

SHWFS 

Shack-Hartmann  WFS 
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3.2.2  ZEMAX  Model. 


ZEMAX  is  a  software  program  used  for  optical  lens  design  and  system  layout. 
Each  optical  surface  is  logged  in  a  spreadsheet  using  characteristics  such  as:  radius 
of  curvature,  propagation  distance,  and  propagation  medium.  A  point  source  (on- 
or  off-axis)  can  then  be  propagated  through  the  optical  layout  using  any  number 
of  geometric  rays.  For  the  purpose  of  this  research,  ZEMAX  was  used  to  model 
the  bench-top  setup  to  determine  the  location  of  the  aperture  stop,  held  stop  and, 
ultimately,  the  magnihcation  of  the  system.  These  details  were  important  to  scale 
the  software-modeled,  elongated  beacon  to  the  bench-top  system. 

The  ZEMAX  model  of  the  system  was  simplihed,  including  only  the  ATS,  the 
4-f  relay,  and  a  single  lenslet  of  the  SHWFS.  Starting  on  the  optical  axis,  a  point 
source  was  hrst  propagated  through  the  model  using  rays.  Incrementally,  the  point 
source  was  placed  further  off  axis,  and  propagated  through  the  model  until  an  optical 
component  began  to  clip  the  rays.  The  hrst  optical  component  to  vignette  the  rays  was 
the  lenslet  while  the  rays  were  only  completely  stopped  at  the  subapertures  detector 
region.  Thus,  these  components  are  the  aperture  stop  and  held  stop,  respectively. 
Figure  35  shows  the  2D  ZEMAX  layout  of  the  lab  bench  where  the  blue  rays  emanate 
from  a  point  source  on  axis  and  the  green  rays  emanate  from  a  point  source  oh  axis. 
Since  the  system  is  very  long  when  shown  with  a  1:1  aspect  ratio,  only  the  ATS  and 
SHWFS  lenslet  are  shown. 

To  hnd  the  theoretical  transverse  magnihcation  My,  rays  from  a  point  source 
separated  a  small  distance  So  oh  the  optical  axis  were  traced  through  the  system. 
These  rays  then  converged  to  a  point  on  the  detector  plane  a  small  distance  oh 
the  optical  axis  S'j.  The  theoretical  transverse  magnihcation  Mt  and  longitudinal 
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Figure  35.  The  ZEMAX  model  of  the  AFIT  lab  bench-top  AO  system.  Only  the  ATS 
and  SHWFS  lenslet  are  displayed  due  to  space  constraints.  The  optical  component 
labels  are  consistent  with  those  in  Table  5.  The  blue  rays  are  from  a  point  source  on 
axis  and  the  green  lines  are  from  a  point  g^urce  off  axis.  The  field  stop  and  aperture 
stop  are  circled  in  red.  REDO  PICTURE!!! 
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Figure  36.  Determining  system  magnification  by  taking  the  ratio  of  image  ray  height 
Si  and  source  ray  height  So  from  the  optical  axis. 

magnification  Ml  was  calculated  by [22]: 


(44) 


Ml  =  M|. 


(45) 


Figure  36  shows  a  simplified  picture  of  how  to  determine  the  theoretical  magnihcation 
of  an  optical  system.  The  magnihcation  parameters  Mt  and  Ml  were  found  to  be 
0.13  and  0.0169,  respectively.  It  should  be  noted  that  magnihcation  determined  by 
ZEMAX  is  only  valid  for  images  formed  in  the  SHWFS.  Lab  experiments,  discussed 
in  the  next  section,  were  used  to  verify  the  ZEMAX  magnihcation. 

3.3  Extended  Source  Model  Progression 

Creating  a  three-dimensional  scattering  volume  on  the  optical  table  was  done  in 
incremental  steps.  First,  diherent  optical  components  were  arranged  to  create  ad¬ 
justable,  extended  2D  sources.  The  extended  sources  were  propagated  through  the 
optical  system  and  captured  by  the  SHWFS  and  tracker  camera.  Proper  measure- 
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ments  were  taken  to  determine  the  actual  magnification  of  the  system  and  verify 
the  ZEMAX  model.  Finally,  the  source  was  tilted  to  produce  source  elongation  and 
imaged  by  the  system. 

Due  to  the  limitations  of  the  bench-top  system  (2mW  maximum  laser  power  and 
fixed  focal  length  of  the  ATS’s  collimating  lens  LI),  determining  the  proper  combi¬ 
nation  of  optical  elements  to  create  an  extended  source  was  an  iterative  process.  The 
main  trade  off  was  making  the  source  diverge  enough  while  still  producing  enough 
irradiance  to  be  detected  by  the  sensors.  To  do  this,  a  method  of  using  an  iris  and 
different  diffusers  is  discussed  in  the  next  section. 

3.3.1  Iris  and  Diffuser. 

To  create  an  adjustable  extended  source  the  1550nm,  2mW  laser  diode  was  first 
attached  to  a  small  collimator  and  carefully  aligned  to  the  ATS’s  optical  axis.  The 
beam  was  then  propagated  through  an  adjustable  iris  placed  at  the  front  focal  point 
of  the  first  lens  LI  in  the  AO  system.  The  light  departing  the  iris  did  not  diverge 
wide  enough  to  fill  the  full  aperture  of  the  SHWFS.  To  correct  this,  a  motor- mounted 
diffuser  was  placed  before  the  iris.  Figure  37  shows  a  picture  of  the  source  setup. 
As  discussed  in  Sec.  2.2,  the  spinning  diffuser  was  necessary  to  reduce  the  spatial 
coherence  in  the  beam.  Figure  38  shows  a  comparison  of  the  image  with  and  without 
the  diffuser  spinning  in  the  tracker  camera.  While  speckle  similarly  affects  both 
sensors,  the  tracker  camera  provides  a  higher-resolution  image  to  depict  the  speckle 
and  lack  thereof. 

Two  collimators  with  different  diameters  (2mm  and  7mm)  and  four  different  dif¬ 
fusers,  each  with  a  different  grid  designation  (P120,  P220,  P600,  and  P1500),  were 
available  for  setting  up  the  source.  The  first  experiments  were  done  to  determine 
which  collimator /diffuser  combination  created  the  most  adjustability  while  still  per- 
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Figure  37.  The  collimated  beam  passes  through  the  iris  to  create  a  diverging,  incoherent 
source.  The  iris  is  then  adjusted  to  the  desired  size. 


Figure  38.  Effects  of  a  diffuser  on  a  collimated,  coherent  source.  Image-plane  speckle 
(left)  is  caused  by  a  stationary  diffuser  and  is  averaged  out  by  a  rotating  diffuser  (right). 
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mitting  enough  irradiance  to  be  captured  by  the  sensors.  Only  one  diffuser  was  used 
at  a  time  due  to  their  high  attenuation.  Two  diffusers,  no  matter  what  grit,  decreased 
the  irradiance  to  undetectable  levels  in  the  sensors.  Section  4.2  displays  the  results  in 
a  table  and  determines  the  best  combination  to  use.  It  should  be  noted  that  the  cor¬ 
rect  combination  was  determined  qualitatively  by  observing  the  irradiance  patterns 
in  the  sensors. 

3.3.2  System  Magnification. 

To  verify  the  magnification  of  the  system,  a  simple  experiment  was  done  using  the 
adjustable  iris  and  spinning  diffuser  setup.  The  iris  (object  So)  was  set  to  0.5mm  and 
its  image  was  captured  by  the  SHWFS.  It  should  be  noted  that  the  actual  diameter 
of  the  source  was  approximately  1mm  since  So  is  defined  as  the  distance  from  the 
optical  axis  to  the  edge  of  the  object.  Using  a  single  subaperture,  the  illuminated 
pixels  across  the  spot  were  counted  and  multiplied  by  the  pixel  pitch  (25/im)  to  get 
the  size  of  the  image  Si.  As  in  the  ZEMAX  calculation,  the  ratio  of  Si  to  So  was 
taken  to  determine  the  transverse  magnification  Mt-  Figure  39  shows  the  SHWFS 
images  formed  by  the  1mm  source.  The  higher-resolution  tracker  could  not  be  used 
for  this  calculation  because  the  two  sensors  may  have  different  fields  of  view,  and  the 
SHWFS  is  the  primary  concern  here.  The  system  magnification  was  determined  to 
be  approximately  0.15  (±0.03  due  to  inaccuracies  with  measuring  the  iris  diameter 
by  hand  and  the  pixel  cut-off  in  the  image).  The  experimental  magnification  matched 
closely  to  the  ZEMAX-determined  magnification  of  0.14. 

3.3.3  Three-Dimensional  Extended  Source. 

Once  the  proper  iris/diffuser  combination  and  magnification  of  the  system  was 
determined,  the  next  step  was  to  produce  extended-source  properties.  It  is  important 


Figure  39.  SHWFS  spots  formed  by  imaging  a  1mm  extended  source.  This  is  a  small 
sample  (16  subapertures)  of  the  entire  SHWFS  array  (81  subapertures).  The  pixels 
encompassed  by  the  white  circle  were  those  used  in  the  magnification  calculation. 

to  mention  that  all  extended  sources  modeled  here  were  not  formed  by  propagating 
light  through  random  media  as  was  done  in  the  software  simulations  due  to  limitations 
with  the  bench-top  setup.  Here,  the  three-dimensional  extended  source  was  modeled 
by  aligning  the  center  of  the  iris  with  the  optical  axis  and  tilting  it  to  different  angles. 
Figure  40  shows  a  picture  of  the  rotated  iris.  Applying  the  system  magnification  de¬ 
termined  in  the  previous  section,  the  beacon  model  must  have  a  diameter  of  1.5mm 
and  a  depth  of  approximately  8mm.  These  dimensions  would  theoretically  create 
similar  SHWFS  image  size  and  depth  as  those  in  the  software  simulations.  Unfortu¬ 
nately,  the  physical  limitations  of  the  iris  and  its  mechanical  mount  prevent  it  from 
achieving  these  dimensions. 

Even  though  this  model  theoretically  cannot  recreate  the  exact  SHWFS  spots  ob¬ 
served  in  the  software  simulations,  it  was  still  used  to  determine  if  spot  elongation 
could  be  simulated  this  way,  on  this  bench-top  setup.  The  iris  was  set  to  90°  (per¬ 
pendicular  to  the  optical  axis),  45°,  and  30°,  and  imaged  by  the  sensors.  The  hole  of 
the  iris  in  each  scenario  was  adjusted  so  that  a  source  width  of  1.5mm  of  light  was 
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Figure  40.  The  collimated  beam  passes  through  the  iris  to  create  a  diverging,  incoherent 
source.  The  iris  is  then  opened  to  adjust  width  and  tilted  to  adjust  beacon  depth. 

always  visible  to  the  system.  This  forced  the  iris  to  be  opened  more  and  more,  which 
was  necessary  to  let  enough  light  through.  Figure  41  shows  a  picture  of  the  three 
different  iris  rotations  and  the  subsequent  increases  in  iris  size.  Each  source  rotation 
was  propagated  through  vacuum  and  turbulence  and  imaged  through  the  sensors  with 
the  AO  loop  closed.  The  turbulence  was  considered  weak  with  a  Rytov  number  of 
0.02  and  an  rg  of  0.77mm  at  the  DM.  With  d/rO  =  0.65,  these  conditions  are  similar 
to  what  was  used  in  the  computer  simulations.  The  data  is  analyzed  in  Sec.  4.2. 
Due  to  manufacturing  challenges  and  time  constraints,  other  methods  of  modeling  a 
three-dimensional  extended  source  were  not  attempted.  Several,  more  sophisticated 
methods  are  suggested  in  the  Ch.  V. 

3.4  Data  Collection  and  Comparison 

For  each  source  experiment,  100  frames  from  the  tracker  camera  and  200  frames 
from  the  SHWFS  were  captured  and  analyzed.  Since  the  diffusers  significantly  re¬ 
duced  the  irradiance  of  the  source,  the  exposure  periods  of  the  sensors  were  relatively 
long.  Table  6  lists  the  sensor  settings  used  to  capture  data  in  all  experiments.  Analysis 
was  done  via  MATLAB  by  reading  in  the  frames,  averaging  the  frames,  and  tempo¬ 
rally  re-sampling  the  data.  Re-sampling  each  set  basically  equalized  the  exposure 


60 


Figure  41.  The  rotation  of  the  iris  to  form  extended  sources  of  different  depths.  The 
depth  of  the  source  increases  as  the  angle  decreases.  In  (a)  the  iris  is  set  perpendicular 
to  the  optical  axis,  giving  the  source  no  depth.  In  (b)  the  iris  is  set  45°  to  the  optical 
axis,  giving  the  source  a  depth  of  1.5mm.  In  (c)  the  iris  is  set  30°  to  the  optical  axis, 
giving  the  source  a  depth  of  2.6mm. 


Table  6.  Sensor  settings  used  to  capture  data  in  all  lab  experiments. 


Capture 

Number  of 

Frame 

Exposure 

Length  [s] 

Frames 

Rate  [Hz] 

Period  [/is] 

Tracker 

0.1979 

102 

515.4 

500 

SHWFS 

0.2062 

200 

970 

150 

period  so  the  data  could  be  accurately  compared. 


3.5  Chapter  Summary 

In  summary,  software  simulations  were  first  used  to  verify  a  documented  method 
for  creating  elongated  beacon  models  in  the  sodium  layer.  These  beacons  were  com¬ 
pared  by  analyzing  the  elongated  spots  they  formed  in  SHWFS  subapertures.  Next, 
a  ZEMAX  model  of  the  AO  system  was  developed  to  determine  the  aperture  stop, 
field  stop,  and  magnification  of  the  system.  On  the  lab  bench-top,  a  simple  experi¬ 
ment  was  done  to  verify  the  magnification  which  ultimately  determined  the  3D  source 
dimensions  to  use  in  the  lab  experiments.  Next,  an  adjustable  iris  and  spinning  dif¬ 
fusers  were  used  on  an  AO  lab  bench-top  to  simulate  perspective  elongation  in  the 


61 


SHWFS.  In  chapter  IV,  simulated  and  experimental  data  are  analyzed  to  determine 
sufficient  methods  of  modeling  perspective  elongation. 
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IV.  Results  and  Analysis 


This  chapter  presents  the  results  and  analysis  for  the  comparison  of  perspective 
elongation  models  developed  in  Ch.  111.  The  results  are  organized  into  three  sections. 
The  first  describes  the  best  way  (according  to  this  study)  to  model  perspective  elon¬ 
gation  with  insight  gained  from  the  computer  simulation  results.  The  second  section 
presents  the  results  from  different  bench-top,  extended-source  experiments.  Finally, 
the  experimental  results  are  compared  with  the  computer  simulation  results. 

4.1  Computer  Simulation  Results 

Perspective  elongation  in  SHWFS  subapertures  was  simulated  by  modeling  an 
elongated  beacon  in  the  sodium  layer  and  propagating  it  through  the  AO  system. 
The  elongated  beacon  was  modeled  19  ways  by  dividing  it  into  different  numbers  of 
evenly-spaced  slices.  This  created  slightly  different  elongated  spot  images  in  each 
SHWFS  subaperture.  Comparing  the  spot  images  by  using  the  techniques  discussed 
in  Sec.  2.9  determined  the  minimum  number  of  slices  needed  to  accurately  model  the 
elongated  beacon. 

4.1.1  RMS  Percent  Difference. 

The  first  method  of  comparing  the  elongated  spots  was  done  by  using  the  RMS 
percent  difference.  Figure  42  shows  the  RMS  percent  difference  of  spot  shapes  formed 
by  adjacent,  increasing  numbers  of  slices  in  both  vacuum  and  turbulence  scenarios, 
using  slices  of  constant  and  column-density-weighted  intensity.  In  all  cases,  the  spot 
shapes  decrease  in  percent  difference  when  formed  by  an  increasing  number  of  slices. 
As  expected,  the  spots  formed  by  column-density-weighted  slices  cause  a  smaller 
initial  difference  due  to  the  relatively  low  intensity  levels  at  the  edges  of  the  sodium 
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Figure  42.  RMS  percent  difference  comparison  of  elongated  spots,  (left)  is  the  percent 
change  in  elongated  spot  shape  formed  through  vacuum  and  (right)  is  the  percent 
change  in  spots  formed  through  turbulence.  Black  represents  sodium  layers  at  constant 
intensity  and  red  represents  layers  with  Gaussian  weighting. 

layer. 

Even  with  a  large  number  of  slices,  the  percent  difference  never  dropped  below  8% 
which  is  too  high  to  dehnitively  determine  how  many  slices  are  adequate.  The  rela¬ 
tively  large  difference  is  probably  because  this  comparison  technique  compares  exact 
data  values  instead  of  how  correlated  the  different  SH  images  are.  Including  another 
two  slices  of  intensity  data  each  simulation  always  increased  the  spot  brightness  in 
the  SHWFS  subapertures.  Another  comparison  technique  to  that  isolates  the  spot 
shapes  with  insensitivity  to  exact  values. 

4.1.2  Correlation  Coefficient. 

A  more  accurate  way  of  comparing  the  change  in  spot  shapes  is  by  calculating  the 
correlation  coefficient  between  corresponding  spot  images.  This  essentially  compares 
the  two-dimensional  flattened  spot  images  formed  by  the  different  number  of  sodium 
layer  slices.  Using  the  39  slice  model  as  the  “continuous”  model,  the  correlation 
coefficient  was  calculated  between  it  and  the  “discrete”  models  (3-37  slices).  Figure 
43  shows  the  plots  of  the  correlation  coefficient  using  slices  of  constant  and  column- 
density-weighted  intensity. 
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Figure  43.  Correlation  coefRcient  comparison  of  elongated  spots,  (left)  are  the  coef¬ 
ficients  for  spots  formed  through  vacuum,  and  (right)  are  the  coefficients  for  spots 
formed  through  turbulence.  Black  represents  sodium  layers  at  constant  intensity  and 
red  represents  layers  with  Gaussian  weighting. 


As  evident  in  the  plot,  the  spot  images  in  both  vacuum  and  turbulence  are  highly 
correlated.  If  we  consider  0.95  to  be  the  correlation  threshold,  then  the  3  slice  elon¬ 
gated  spot  model  is  sufficient  for  the  vacuum  case.  For  the  turbulence  case,  the  9  slice 
elongated  spot  model  is  sufficient.  As  a  reminder,  the  elongated  spots  were  formed 
using  propagations  near  zenith.  As  the  angle  from  zenith  increases,  the  beacon  elon¬ 
gates  more.  In  these  cases,  more  slices  may  be  necessary  to  sufficiently  model  the 
elongated  beacon. 


4.2  Experimental  Results 

The  experimental  results  and  analysis  are  presented  here.  The  proper  diffuser 
and  iris  combination  for  our  AO  bench-top  setup  was  determined  first.  With  the 
proper  combination,  the  iris  was  adjusted  and  rotated  to  create  elongation  in  SHWFS 
subapertures.  The  subsections  below  provide  all  relevant  details  and  imagery. 
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Table  7.  The  divergence  and  irradiance  of  different  diffusers  applied  to  two  collimated 
2mW  beams. 


diffuser  grit  classification 

120  220  600 

1500 

beam  diameter 

mm 

div. 

irrad.  div.  irrad.  div.  irrad. 

div. 

irrad. 

2 

good 

good  good  good  poor  good 

poor 

good 

7 

good 

poor  good  poor  poor  good 

poor 

good 

4.2.1  Iris  and  Diffuser. 

Pairing  each  diffuser  with  the  adjustable  iris,  the  divergence  of  the  source  and 
irradiance  was  observed  by  the  tracker  camera  and  SHWFS.  Table  7  shows  the  per¬ 
formance  of  each  diffuser  using  two  different  collimated  beams.  The  only  combination 
that  produced  a  wide  enough  divergence  angle  and  enough  irradiance  was  the  2mm 
collimator  paired  with  the  coarsest  diffusers  (120  and  220).  The  7mm  beam  would 
have  provided  the  most  flexibility  in  the  source,  but  the  power  was  too  weak  once  it 
was  spread  that  much.  All  other  combinations  either  did  not  produce  enough  diver¬ 
gence  or  enough  irradiance.  All  observations  were  qualitative.  While  this  experiment 
did  not  provide  any  direct  contribution  to  producing  perspective  elongation,  it  was  a 
key  step  in  the  process  and  gives  the  boundaries  for  modeling  an  extended  source  on 
this  lab  bench-top.  Ultimately  the  2mm  collimator  with  the  120  diffuser  was  chosen. 


4.2.2  Three-Dimensional  Extended  Source. 

Using  the  illuminated  adjustable  iris  and  diffuser  as  the  source,  the  iris  was  ro¬ 
tated  to  three  different  positions  to  create  three  different  extended  sources.  The  first 
position  was  with  the  iris  at  90°  (flat),  the  second  was  with  the  iris  at  45°,  and  the 
third  was  with  the  iris  at  30°.  The  three  source  positions/shapes  were  captured  by  the 
tracker  camera  and  the  SHWFS  after  propagating  through  vacuum  and  turbulence. 
Figure  44  shows  the  sensor  images  through  vacuum,  and  Fig.  45  shows  the  sensor 
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images  through  weak  turbulence.  As  the  angle  of  the  iris  decreased,  the  images  in  the 
sensors  became  more  elliptical  as  was  expected.  In  Figs.  44  and  45  the  tracker  image 
and  corresponding  SHWFS  images  are  paired.  The  top  pair  of  images  (a)  is  with  the 
iris  at  90°,  the  middle  pair  (b)  is  with  the  iris  at  45°,  and  the  bottom  pair  (c)  is  with 
the  iris  at  30°.  In  the  bottom  pair  of  images,  a  larger  region  of  interest  (ROI)  was 
needed  to  display  the  elongated  beacon  on  the  tracker.  This  makes  the  image  seem 
smaller  than  the  other  two,  but  careful  inspection  of  the  dimensions  show  it  is  on  the 
same  order. 

The  vacuum  and  turbulence  images  are  nearly  indistinguishable.  This  is  mainly 
because  a  small  d/ro  value  (approximately  0.649)  was  used  in  the  turbulence  scenario. 
In  this  scenario,  the  spots  did  not  show  additional  spreading,  but  did  display  a  slight 
tilt.  However,  since  all  images  were  captured  in  the  closed-loop  regime,  the  tilt  is  not 
apparent  in  Fig.  45.  A  larger  d/r^  would  most  likely  spread  the  spots  much  more 
and  introduce  more  tilt,  but  this  was  not  not  tested  here  and  is  outside  the  scope  of 
this  research.  Additionally,  while  not  part  of  the  investigation,  it  was  observed  that 
the  closed-loop  Strehl  ratio  decreased  with  increasing  spot  size,  as  expected. 

The  SHWFS  spots  formed  by  the  extended  source  model  cannot  be  quantitatively 
verihed  with  data  from  the  software  simulations.  This  is  because  the  source  dimen¬ 
sions  to  properly  model  the  sodium  beacon  could  not  be  achieved  by  just  simply 
rotating  the  iris.  During  the  course  of  the  experiment,  it  was  discovered  that  match¬ 
ing  the  depth  of  the  beacon  was  outside  the  physical  limitations  of  the  simple  iris.  The 
exact  iris  dimensions  would  have  been  1.5mm  by  20mm,  forcing  it  to  be  at  a  4°  angle 
with  the  optical  axis.  Ultimately,  the  narrow  angle  decreased  the  throughput  of  the 
beam  so  much  that  the  source  could  not  be  imaged.  More  robust  source  designs  must 
be  developed  in  order  to  meet  the  required  dimensions,  and  are  outlined  in  Ch.  5. 
However,  by  doing  the  in-depth  magnification  experiments  and  calculations,  the  spots 
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are  beginning  to  approach  the  sizes  expected  to  be  seen  at  SOR.  Unfortunately,  even 
a  direct  qualitative  comparison  cannot  be  accomplished  for  two  reasons.  First,  SOR 
data  was  not  releasable  at  the  time  of  this  research.  Second,  and  more  importantly, 
SOR  uses  a  quad-cell  based  SHWFS  and  tracker.  Essentially,  each  subaperture  only 
consists  of  four  pixels  so  information  about  spot  shape  is  not  captured,  while  the  AO 
system  in  our  lab  has  high-resolution  subapertures  and  a  full  array  tracker  which  has 
enough  resolution  to  characterize  the  spot  shape. 

4.3  Chapter  Summary 

This  chapter  presented  the  results  from  both  the  software  simulation  of  extended 
beacons  and  the  lab  experiments  using  extended  sources.  In  all  cases,  the  extended 
beacon  models  were  evaluated  by  analyzing  the  spot  shapes  they  formed  in  SHWFS 
subapertures.  The  software  model  showed  that  dividing  the  elongated  beacon  into  3 
slices  (for  propagation  through  vacuum)  and  9  slices  (for  propagation  through  turbu¬ 
lence)  adequately  modeled  the  elongated  beacon  in  the  sodium  layer.  This  Ending  is 
consistent  with  claims  made  in  current  literature  [5]. 

The  bench-top  extended  source  demonstrated  preliminary  steps  toward  modeling 
perspective  elongation  in  a  lab  setting.  By  rotating  and  adjusting  an  iris  illuminated 
incoherently  by  laser  light  passing  through  a  moving  diEuser,  various  levels  of  elon¬ 
gation  were  observed  in  SHWFS  subapertures.  It  should  be  stressed  that  it  is  not 
clear  if  all  properties  of  perspective  elongation  were  created  by  the  model,  but  the 
adjustable  model  is  useful  for  many  areas  of  future  research  as  discussed  in  Ch.  5. 
While  much  more  work  needs  to  be  done  to  accurately  model  perspective  elongation 
on  a  lab  bench-top,  the  experiments  here  will  help  pave  the  way  by  providing  the  first 
documented  tests  of  extended  sources  on  bench-top  AO  systems. 
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Figure  44.  Extended  source  and  resultant  elongated  SHWFS  spots  through  vacuum. 
Starting  with  the  top  pair  of  images,  (a)  is  with  the  iris  set  perpendicular  to  the  optical 
axis  (90°),  (b)  is  with  the  iris  set  to  45°,  and  (c)  is  with  the  iris  set  to  30°. 
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Figure  45.  Extended  source  and  resultant  elongated  SHWFS  spots  through  turbulence. 
Starting  with  the  top  pair  of  images,  (a)  is  with  the  iris  set  perpendicular  to  the  optical 
axis  (90°),  (b)  is  with  the  iris  set  to  45°,  and  (c)  is  with  the  iris  set  to  30°.  All  images 
were  captured  in  the  closed-loop  regime. 
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V.  Conclusions  and  Recommendations 


This  chapter  summarizes  all  pieces  of  this  research.  Key  results  from  data  analysis 
in  Ch.  IV  are  listed  for  easy  reference.  Additionally,  recommendations  for  future  work 
on  perspective  elongation  are  suggested. 

5.1  Summary 

As  discussed  in  Ch.  I,  the  objectives  of  this  research  were  to  validate  a  method  for 
modeling  an  elongated  beacon  with  wave  propagation  software,  perform  lab  experi¬ 
ments  to  simulate  perspective  elongation  in  SHWFS  subapertures,  and  deliver  work¬ 
ing,  scalable  software  and  lab  models  for  simulating  perspective  elongation.  Chapter 
II  provided  all  the  necessary  background  information  used  throughout  the  research. 
In  Ch.  Ill,  software  and  lab  environments  were  developed  to  meet  Ch.  I  objectives. 
Using  wave  propagation  software  in  MATLAB,  an  atmospheric  turbulence  model  was 
developed.  Gaussian  beams  were  propagated  through  the  model  to  the  sodium  layer 
to  create  an  elongated  beacon.  Dividing  the  sodium  layer  into  different  slices  created 
slightly  different  elongated  beacon  models.  The  models  were  analyzed  by  imaging 
them  with  a  SHWFS  and  comparing  the  resultant  elongated  spots.  The  lab  portion 
involved  trying  different  methods  of  modeling  an  elongated  scattering  source  to  get 
elongated  SHWFS  spots.  Two-dimensional  sources  were  used  to  first  determine  the 
source  width  and  a  three-dimensional  source  was  used  to  determine  the  source  depth. 
Chapter  IV  analyzes  the  data  from  software  and  lab  experiments.  Data  from  the 
software  model  was  used  to  verify  the  lab  model.  The  next  section  summarizes  the 
novel  contributions  from  the  software  and  lab  experiments. 
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5.2  Novel  Contributions 


The  novel  contributions  obtained  in  this  study  are: 

•  Software  model  that  recreates  perspective  elongation.  The  software  program 
effectively  simulated  the  propagation  of  a  sodium  laser  source  to  create  an  elon¬ 
gated  beacon.  The  beacon  was  then  imaged  by  a  SHWFS  to  simulate  elongated 
spots.  Several  parameters  were  adjustable  such  as:  separation  distance  between 
the  primary  aperture  and  laser  launching  source,  all  propagation  distances,  tur¬ 
bulence  strength,  and  SHWFS  dimensions  (number  of  lenslets,  focal  length, 
lenslet  dimensions). 

•  Number  of  slices  needed  to  sufficiently  model  and  elongated  sodium  beacon.  The 
number  of  slices  to  sufficiently  model  an  elongate  beacon  was  determined  by  an¬ 
alyzing  the  resulting  SHWFS  spots.  In  vacuum  propagation,  it  was  determined 
that  only  3  slices  were  needed  to  accurately  model  the  elongated  beacon.  In 
turbulence  propagation,  it  was  determined  that  at  least  9  slices  were  needed. 

•  Bench-top  extended  source  model.  The  bench-top  source  was  created  by  propa¬ 
gating  laser  light  through  a  moving  diffuser  and  an  adjustable  iris.  The  moving 
diffuser  was  used  to  make  the  light  spatially  incoherent  for  imaging  purposes. 
The  iris  was  was  adjusted  (opened  and  turned)  to  create  the  desired  width  and 
depth  of  the  extended  source. 

•  Added  capabilities  to  work  with  extended  beacons.  While  theorized  about,  the 
effects  of  extended  beacons  on  WFS’s  have  been  widely  untested.  With  the  ba¬ 
sic  source  setup  used  in  this  research,  wavefront  sensing  with  correlation-based 
slope  measurements  can  be  analyzed  for  different  WFS  such  as  the  SHWFS  and 
the  Self- Referencing  Interferometer  (SRI).  Additionally,  the  extended  sources 
can  be  used  to  validate  the  effects  of  deep  turbulence  on  AO  systems [23]  [24]  [25]. 
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Deep  turbulence  is  essentially  an  extended  source  through  strong  turbulence. 
Strong  turbulence  has  been  well  tested,  but  not  when  combined  with  an  ex¬ 
tended  source.  Thus,  the  model  developed  here  adds  the  hnal  piece  to  test  deep 
turbulence  theory. 

5.3  Recommendations 

This  section  provides  recommendations  for  future  work  related  to  this  research. 

•  Elongated  beacon  software  model.  The  elongated  spots  used  in  the  software 
simulation  were  modeled  by  dividing  the  sodium  layer  into  slices.  In  each  case, 
the  slices  were  always  equally  spaced  between  80  and  100km  and  centered  at 
90km.  A  better  method  may  be  to  £x  a  slice  at  90km  and,  keeping  the  distance 
constant,  increase  the  number  of  slices  equally  above  and  below  the  center  slice. 
Figure  46  shows  a  picture  of  this.  This  method  may  produce  different  results 
for  the  column-density-weighted  scenario  since  the  sodium  atom  concentration 
is  greatest  around  90km.  For  the  constant  weighted  scenario,  the  results  would 
be  close  if  not  exactly  the  same. 

•  Telescope  and  atmospheric  parameters  in  software  simulations.  All  software 
simulations  were  done  modeling  a  3.5m  telescope  pointed  at  zenith  in  relatively 
weak,  static  turbulence  conditions.  Changing  any  and/or  all  of  these  variable 
will  dehnitely  have  an  impact  on  the  results.  Assumedly,  harsher  conditions 
(larger  telescope,  angles  from  zenith,  and  stronger,  dynamic  turbulence)  will 
need  more  slices  to  adequately  model  an  elongated  beacon. 

•  Extended  source  lab  models.  In  this  study,  the  bench-top  extended  source  model 
was  created  by  illuminating  a  rotated  iris  with  incoherent  light.  While  this 
created  perspective  elongation  in  the  SHWFS  subapertures,  it  had  limits  with 
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Figure  46.  Two  different  ways  of  dividing  the  sodium  layer  to  model  elongated  spots. 
The  top  figure  shows  the  method  used  in  this  report  and  the  bottom  figure  shows  a 
possible  future  setup.  A  3,  5,  and  7  slice  model  is  given  for  both  setups. 
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changing  the  dimensions  of  the  sonrce.  Three  other  ways  to  model  a  three- 
dimensional  extended  beacon  are  presented  here. 

The  hrst  model  called  the  ’’rice  grain”  model  involves  nsing  a  highly  reflective 
object  ronghly  the  size  of  a  rice  grain.  Monnted  on  the  rotating  shaft  of  a 
motor  and  aligned  with  the  optical  axis,  the  grain  wonld  be  illnminated  by 
a  coherent  laser  beam.  Using  the  motor  to  rotate  the  grain  would  create  an 
incoherent  source  to  be  imaged  by  the  system.  The  rice  grain  could  be  rotated 
to  any  angle  to  create  an  extended  three-dimensional  source  that  matches  the 
demagnihed  dimensions  of  an  elongated  sodium  beacon.  Figure  47(a)  shows  a 
picture  of  the  ’’rice  grain”  model. 

The  second  model  called  the  ’’reflector  tube”  model  is  very  similar  to  the  ’’rice 
grain”  model  in  that  is  uses  a  devices  mounted  on  a  motor  and  aligned  on  the 
optical  axis  of  the  AO  system.  This  source  is  made  by  injecting  a  clear  resin 
mixed  with  highly  reflective  powder  into  a  small  glass  tube.  The  device  is  then 
mounted  on  a  motor  shaft  and  illuminated  by  a  laser.  The  tube  could  be  rotated 
to  any  angle  to  create  an  extended  three-dimensional  source  that  matches  the 
demagnihed  dimensions  of  an  elongated  sodium  beacon.  47(b)  shows  a  picture 
of  the  ’’rehector  tube”  model. 

The  third  model  called  the  ’’LED  bar”  model  involves  mounting  several  tiny 
light  emitting  diodes  (LED)s  on  an  adjustable,  mechanical  mount.  Each  diode 
would  face  the  same  direction,  all  pointed  along  the  optical  axis.  The  center 
of  the  bar  would  be  positioned  on  the  optical  axis  at  the  focus  of  the  entrance 
pupil,  and  could  be  rotated  to  any  position.  Also,  the  width  of  the  bar  could 
be  adjusted  by  expanding  or  contracting  the  device.  Figure  48  shows  a  picture 
of  the  ’’LED  bar”  model. 
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Figure  47.  Two  devices  to  possibly  model  and  extended  source.  The  rice  grain  (a)  and 
the  reflective  tube  (b)  are  both  aligned  with  the  optical  axis  and  mounted  on  a  motor. 
A  laser  is  reflected  off  the  spinning  sources  which  could  be  rotated  to  any  desired  angle. 


Figure  48.  A  deformable  LED  bar  to  possibly  model  an  extended  source.  The  bar 
could  be  expanded  or  contracted  and  rotated  along  any  axis. 
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5.4  Chapter  Summary 


This  chapter  summarized  the  key  results  and  presented  some  future  research  rec¬ 
ommendations.  The  software  model  verified  a  documented  beacon  model  to  recreate 
perspective  elongation,  and  the  bench-top  source  model  was  a  solid  first  step  in  simu¬ 
lating  perspective  elongation  in  the  lab.  Future  recommendations  included  a  different 
software  model  to  simulation  perspective  elongation  and  a  few  different  bench-top 
sources  to  better  model  an  extended,  three-dimensional  source. 
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Appendix  A.  MATLAB  Code 


Listing  A.l.  LGS_Source.m 

1  %  Grid  sampling  parameters 

2  N  =  1024;%  512  number  of  grid  points 

3  An  =  0.01  ;  %  0.02  grid  spacing  in  sodium  layer  [m] 

4  Al  =  0.004;  %  0.004  grid  spacing  in  source  plane  [m] 

5  L  =  N*An;  %  actual  length  of  grid 

6  wvl  =  589. 2e- 9;  %  wavelength  of  sodium  guidestar 

7  k  =  2*  pi /wvl;  %  wave  number 

8  z  =  linspace  (85e3  ,  95e3  ,  nlayers);  %  propagation  distance 

9  %  source  parameters 

10  W_0  =  0.08;  %  diameter 

11  F_0  =  90000;  %  curvature 

12  Theta.O  =  l-z/F_0;  %  focus 

13  Lambda_0  =  2  .  *  z  /  ( k^^W.O  "  2 ) ;  %  diffraction 

14  %  receiver  parameters 

15  W  =  W-0. *  sqrt  (  Thet  a_0 .  "  2  +  Lambda_0.  "  2 ) ;  %  diameter 

16  F  =  F_0*(Theta_0.  ''2+Lambda_0.  ''2).*(Theta_0-l)  ... 

17  .  /  (  Theta_0 . ''2+Lambda_0. "  2  -  Theta.O  ) ;  %  curvature 

18  Theta  =  Theta_0. /(  Theta_0.  "2  +  Lambda_0.  "  2  ) ;  %  focus 

19  Lambda  =  Lambda-0. /(  Theta_0.  "2  +  Lambda-0.  "  2 ) ;  %  diffraction 

20  %  cartesian  and  polar  coords  (source  plane) 

21  [xl  yl]  =  meshgrid  ( ( -N/2  :  N/2-1)*a1);  %  grid 

22  [Thl  rl]  =  c ar t 2 pol ( xl , y 1 ) ;  %  convert  to  polar  coords 

23  %  Gaussian  beam  field  and  irradiance  at  source 

24  U_gaus_s  =  exp  (  -  (  r  1 .  "  2 )  .  /  (W_0 "  2 ) )  .  *  exp  ( -  i  *  k .  *  (  r  1 .  "  2)  .  /  ( 2*  F_0  ) )  ; 

25  I_gaus_s  =  abs(U_gaus_s)."2; 


Listing  A. 2.  LGS  J*rop Vac.ni 

1  %  load  vacuum  ground  to  sky  (g2s)  model 

2  load  '  Turb.models  ' 

3  z  =  linspace  (80e3  ,  100e3  ,  7);  %  propagation  distance  (7  slices) 

4  zm  =  mean(z);  %  middle  Na  layer 

5  LGS-Source  %  load  Gaussian  beam  source  parameters 
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6 

7  %  cartesian  and  polar  coords  (receiver  plane) 

8  [x2  y2]  =  meshgrid((-N/2  :  N/2-l)*An);  %  grid 

9  [Th2  r2]  =  cart2pol {x2 , y2 ) ;  %  convert  to  polar  coords 
10 

11  E_gaus  =  emfield  ( U_gaus_s  ,  wvl  ,  N*a1);  %  Make  field  for  WaveProp 

12  Vac _g2 s  .  sc aling^An/A  1 ;  %  Proper  vac  scaling  factor  for  WaveProp 

13 

14  %  Propagate  Gaussian  field  through  vac/turb  model 

15  Ev_gaus_g2s  =  Vac_g2s  (  E.gaus  )  ; 

16 

17  z2  =  z-50e3;  %  distance  from  50km  to  7  layers  in  mesosphere 

18  layers  =  7;  %  vac  propagation  to  7  slices  in  sodium  layer 

19  emarray_Na_v7  =  zeros  (  [N  N  layers  ]  )  ; 

20  for  idx  =  l:  layers  ;  %  index  for  propagation  distance 

21  V2  =  vacuum  ( z2  ( idx  ))  ;  %  propagate  to  the  sodium  layer 

22  Ev_gaus_g2s_2  =  V2 (  Ev_gaus_g2s  )  ;  %  WaveProp  field  thru  vacuum 

23  Iv-gaus  =  intensity  ( Ev_gaus_g2s_2  ) ;  %  WaveProp  intensity  thru  vacuum 

24  %  Theoretical  Gaussian  beam  field  at  receiver  (from  A&P) 

25  U.gaus.r  =  l./(Theta_0(idx)+i.  *Lambda_0  ( idx  ) )  .  *  exp  (i.*k.*z(idx))  ... 

26  . *  exp (-(r2."2)./  (W(  idx)  .  "2))  .*  exp (-i*k.*(r2.  "2)  ./(2.*F(idx))); 

27  %  Theoretical  Gaussian  beam  irradiance  at  receiver  (from  A&P) 

28  I_gaus_r2  =  abs  (  U -gaus  _r  )  .  ''2\ 

29  end 


Listing  A. 3.  MeanJrradiance.m 

2  WaveProp  models  of  Gaussian  beam  propagation  thru  Turbulence  %%% 

4  load  '  Turb.models  '  %  load  turbulence  model 

5  Turb_g2s.proprange  =  90e3;  %  set  propagation  distance  in  turb  model 

6  Tur  b  _g2s  .  alt  i  t  ude  2  =  90e3;  %  set  altitude  in  turb  model 

7 

8  z  =  90e3;  %  constant  propagation  distance 

9  LGS-Source  %  load  Gaussian  beam  source  parameters 
10 

11  %  cartesian  and  polar  coords  (receiver  plane) 


79 


12  [x2  y2]  =  meshgrid  ( ( -N/2  :  N/2-l)*An);  %  grid 

13  [Th2  r2]  =  cart2pol (x2 , y2 ) ; 

14 

15  E_gaus  =  emfield  ( U_gaus_s  ,  wvl  ,  N*Al);  %  Make  field  for  WaveProp 

16  Turb_g2s.scaling^\n/Al ;  %  Proper  scaling  factor  for  WaveProp 

17 

18  1 1  _g  aus  _g2  s  _2  =  zeros(N);  %  initialize  variable  in  loop 

19  props  =  10000;  %  number  of  propagations 

20  for  idx  =  1:  props 

21  Turb_g2s.seed  =  r andi  ( 1 0 0 00 ) ;  %  change  seed  every  prop 

22  Et_gaus_g2s  =  Turb_g2s  ( E_gaus  )  ;  %  propagation  thru  turb  model 

23  It_gaus_g2s  =  intensity  (  Et_gaus_g2s  )  ;  %  intensity 

24  It_gaus_g2s_2  =  1 1  _gaus  _g  2  s  _2  +  It_gaus_g2s;  %  sum  each  irrad  pattern 

25  end 

26  %  WaveProp  model  of  turb  degraded  irradiance  (averaged) 

27  I.gaus.turb  =  It_gaus_g2s_2 /props ; 

28 

29 

30  Analytical  Turbulence  Degraded  Irradiance  %%% 

31 

32  h  =  0:0.05:z;  %  Step  size  for  integral 

33  zen  =  0*pi/180;  %  zenith  angle  [rad] 

34  %  HV-57  Cn2  Model 

35  v  =  21;  %  velocity  for  HV57  [m/ s  ] 

36  A=1.7e-14;  %  Coefficient  for  HV57 

37  Cn2  =  0. 00594 *(v/27)'2*(h*le-5). MO . *exp ( -h/1000) . . . 

38  +2 . 7e  -  16*  exp  (  -  h/1500)+A*exp  (  -  h/ 1 00 ) ; 

39  %coherence  length  (spherical  wave)  to  test 

40  rO-down  =  (0.423*k"2*trapz(Cn2.*((z-h)./z)."(5/3))  ... 

41  *  ( h  ( 2 )  -  h  ( 1 ) )  *  sec  ( zen  ) )  "  (  -  3/ 5) 

42  %  calculate  variance  factor 

43  sigma2r  =  8.3. *((  Lambda  ( idx2)*z)/(k*  rO.down  ^  2 ) )  ... 

44  .  "(5/6)  .  *((  r2.  M)  .  /(W(idx2  )  "2))  ; 

45  %  calculate  T 

46  T=9.98*((  Lambda  ( idx2)*z)/(k*  rO.down  "2))  .  "(5/6); 

47  %  analytical  model  of  turb  degraded  irradiance  (A&P  page  189,  Eqs.  45  &  46) 

48  WT.T  =  W(idx2)*sqrt(l+T); 

49  I_gaus_turb_an  =  ( ( W-0 "  2 )  /  (WILT "  2 ) )  .  *  exp  ( -  ( 2*  r  2  .  "  2)  .  /  (WILT "  2  ) ) ; 
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Listing  A. 4.  FulLPropagation.m 


1 

2  Wo%  Simulating  Perspective  Elongation  in  SHWFS  Subapertures  %%% 

3  1  .  Beam  propagation  to  sodium  layer  %%% 

4  2.  Beacon  propagation  to  telescope  from  sodium  layer  %%% 

5  3.  Propagtions  done  through  both  vacuum  and  turbulence  %%% 

7 

8  load  '  Turb.models  '  %load  vac  and  turb  g2s/s2g  models 

9 

10  vbase  =  ' Ivspots_nlayers_ ' ;  %  base  variable  name  for  uncombined  vac  spots 

11  tbase  =  '  Itspots_nlayers_  '  ;  %  base  varialbe  name  for  uncombined  turb  spots 

12  nlayers  =  [3,5,7,9,11,13,15,17,19,21,23,25,27,29,... 

13  31,33,35,37,39];  %  Na  slices 

14  Ispots_v  =  zeros([560  560  length  ( nlayers  )])  ;  %  initialize  spot  turb  vars 

15  Ispots.t  =  zeros([560  560  length  ( nlayers  )])  ;  %  initialize  spot  vac  vars 

16  for  idxl  =  1 :  length  (  nlayers  )  ;  close  all;  %  choose,  loop  over  slice  models 

17  z  =  linspace  (80e3  ,  100e3  ,  nlayers ( idxl )) ;  %  propagation  distance 

18  zm  =  mean(z);  %  middle  Na  layer 

19  LGS-Source  %  calls  gaussian  beam  source  parameters 

20 

21  LGS-sep  =  4;  %  LGS  separation  from  center  of  tele  ap  [m] 

22  tseed  =  666;  %  constant  turbulence  seed  for  up  Si  down  props 

23 

24  [xn  yn]  =  meshgrid  ( ( -N/2  ;  N/2-l)*An);  %  grid  in  receiver  planes 

25  [Thn  rn]  =  cart2pol (xn , yn ) ; 

26 

27 

28  Propagate  from  ground  to  Na  slices  %%%%% 

29 

30  E.gaus  =  emfield  ( U_gaus_s  ,  wvl  ,  N*a1);  %  make  gaussian  field  for  WP 

31 

32  V ac _g2 s . scali ng  =  An/Al;  %  proper  scaling  for  WP 

33  Turb_g2s. scaling  =  An/Al;  %  Proper  scaling  for  WP 

34  Turb_g2s.seed  =  tseed  ;  %  set  the  turb  seed 

35 

36  Ev_gaus_g2s  =  Vac_g2s  ( E.gaus  ) ;  %  prop  gauss  fid  50km  thru  vac  model 

37  Et_gaus_g2s  =  Turb_g2s  ( E.gaus  )  ;  %  prop  gauss  fid  50km  thru  turb  model 

38 

39  z.up  =  z-50e3;  %  distance  from  50km  to  incremental  NA  slices 
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40  Iarray_Na_v  =  zeros([NN  nlayers(idxl)]);  %  initialize  vac  up  prop 

41  Iarray_Na_t  =  zeros([NN  nlayers(idxl)]);  %  initialize  turb  up  prop 

42  matlabpool  ;  %  initialize  the  matlab  pool 

43  for  idx  =  1 :  nlayers  (  idxl  )  ;  %  loop  thru  each  slice; 

44  V_up  =  vacuum(  z_up  ( idx  ) )  ;  %  set  vac  distance  from  50km  to  slices 

45  Ev_gaus_g2s_2  =  V_up  (  Ev_gaus_g2s  )  ;  %  prop  gauss  beam  thru  vac 

46  Et_gaus_g2s_2  =  V_up  (  Et  _gaus _g2s  )  ;  %  prop  gauss  beam  thru  turb 

47  Iv_gaus  =  intensity  ( Ev_gaus_g2s_2  )  ;  %  find  vac  irradiance 

48  It.gaus  =  intensity  (  Et_gaus_g2s_2  )  ;  %  find  turb  irradiance 

49 

50  %  store  layers  (irradiance  ONLY)  in  frames 

51  Iarray_Na_v(:,;,idx)  =  [Iv.gaus.data]; 

52  Iarray_Na_t(:,:,idx)  =  [It.gaus. data]; 

53  end 

54 

56  Propagate  from  Na  slices  to  telescope  aperture  %%%%% 

58  A50km  =  0.0036;  %  grid  spacing  at  50km  [m] 

59  z_dn  =  z-50e3;  %  distance  from  levels  in  Na  layer  to  50km 

60  [x50km  y50km]  =  meshgrid  ( ( -N/2  :  N/2 -  1 ) * A50km ) ;  %  grid  at  50km 

61 

62  %  super  gaussian  (for  filtering) 

63  sigma  =  .15*L;  n  =  16; 

64  SG_x  =  exp  (  -  ( ( x50km )  /  sigma )  .  "  n  )  ; 

65  SG-y  =  exp  (-( (y50km)  /  sigma ).  "n  )  ; 

66  SG-S  =  SG-x.*SG-y;  %  3D  square  ’’gaussian” 

67 

68  Vac_s2g . scaling  =  An/A50km;  %  WP  scaling  at  50km 

69  Turb_s2g. scaling  =  An/A50km;  %  WP  scaling  at  50km 

70  Turb_s2g.seed  =  tseed  ;  %  set  the  turb  seed  (same  as  above) 

71  emarray _ap_v  =zeros([NN  nlayers(idxl)]);  %  initialize  vac  down  prop 

72  emarray_ap_t  =  zeros  ([NN  nlayers(idxl)]);  %  initialize  turb  down  prop 

73  for  idx  =  1 :  nlayers  (  idxl  )  ;  %  loop  thru  each  slice; 

74  xc  =  (  z  ( idx  ) -zm)  /  (zm/LGS_sep  ) ;  %  sph  wave  x  shift  from  op  axis 

75  yc  =  -  (  z  ( idx  ) -zm)  /  (zm/LGS_sep  ) ;  %  sph  wave  y  shift  from  op  axis 

76 

77  U-sph_an  =  ( 1  . /z_dn  ( idx  ) )  ...  %  spherical  wave  equation 

78  .  *  exp  (li*(k./(2.*z_dn(idx)))  .4;(  (x50km-  xc  )  .  "  2-|-(y50km-  yc  )  .  "  2  ) ) ; 

79  E_sph  =  emfield  ( U_sph_an  ,  wvl  ,  N*A50km)  ;  %  make  EM  field  for  WP 
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E.sph  =  E_sph.*SG_s;  %  filter  field  with  super  gaussian 


81 

82  Ev_sph_s2g  =  Vac_s2g  ( E_sph  )  ;  %  prop  sph  wave  through  vac 

83  Et_sph_s2g  =  Turb_s2g  ( E_sph  )  ;  %  prop  sph  wave  through  turb 

84 

85  %  store  fields  at  telescope  aperture  in  frames 

86  emarray_ap_v(:,:,idx)  =  [Ev_sph_s2g.data]; 

87  emarray  _ap_t  (:,:,idx)  =  [Et_sph_s2g.data]; 

88  end 

89 

90 

91  9^0%%%  Make  lenslet  array  and  subaps  %%%%% 

92 

93  nsubs  =  40;  %  number  of  subaps  per  side  (based  on  rO) 

94  d  =  D/nsubs;  %  diameter  of  each  lenslet  ,  (based  on  rO) 

95  ripix  =  16;  %  number  of  pixels  per  subap 

96  dpix  =  d/npix;  %  diameter  of  each  pixel  in  subap 

97  foe  =  0 . 1 2  5  *  ( d*  nsubs  *  dpix  )  /  ( 2*  wvl ) ;  %  focal  length  of  lenslets 

98  Llets  =  lensletelement(d,foc);  %  make  lenslets  in  WP 

99  V.Llets  =  vacuum(foc);  %  prop  distance  from  lenlets  to  detector 
100 

Loi 

102  Make  elongated  spots  %%%%% 

L03  WMWWWSWOTWWWWWwm 

104  %  cartesian  coords  for  for  tele  ap 

105  Nt  =  200;  %  grid  pts  across  tele  ap 

106  [xt  yt  ]  =  meshgrid  ( ( -Nt/2  :  Nt/2-l)*An);  %  grid 

107 

L08  %  field  propagated  from  middle  Na  layer  (convert  back  to  emfield) 

109  %  for  collimating  the  fields 

110  mid  =  ceil((length(  emarray  _ap_v  (l,l,:)))/2);  %  find  middle  frame 

111  Ev.mid  =  emfield  (  emarray _ap_v  ( :  , :  ,  mid  )  ,  wvl  ,  L  )  ;  %  make  into  field 

112  Ev.mid  =  cookiecut  ( Ev.mid  ,D) ;  %  cut  out  a  window  size  of  tele  ap  (2m) 

113  Ev.mid.ph  =  phasor  ( Ev.mid  )  ;  %  keeps  phase  only 

114 

115  %  grid  spacing  to  increase  resoltion  of  field  at  tele  ap 

116  An2  =  An  /  2 ; 

117  Nt2  =  Nt*2; 

118  [xt2  yt2]  =  meshgrid  ((- Nt2  /  2 :  Nt2/2  -  1 )  *An2  )  ; 

119 


83 


L20  %  grid  spacing  for  demagnified  field  (field  in  Na  layer) 

121  Aa  =  0.05 /200; 

122  [xa  ya]  =  meshgrid  ( (  -  200/ 2 : 200/2  -  1 )  *  Aa  ) ; 

123  %  grid  spacing  to  decrease  resolution  of  field  in  Na  layer  (for  conv) 

124  Ab  =  0.05/14; 

125  [xb  yb]  =  meshgrid  ( (  -  1 4/ 2 : 1 4/ 2  -  1 )  *  Ab  )  ; 

126 

127  Iv_final  =  zeros(560);  %  initialize  final  irradiance  (vac) 

128  It.final  =  zeros(560);  %  initialize  final  irradiance  (turb) 

129  for  idx  =  1 :  nlayers  (  idxl  )  ;  %  loop  thru  each  slice; 

130  manipulate  fields  at  tele  ap  (sky  to  ground  prop) 

131  Ev  =  emfield(  emarray_ap_v  ( :  , :  ,  idx  )  ,  wvl  ,  L);  %  convert  back  to  fid 

132  Et  =  emfield(emarray_ap_t(:,:,idx),  wvl  ,  L);  %  convert  back  to  fid 

133  %  cut  out  2nix2m  window 

134  window_Ev  =  cookiecut  (Ev,D)  ;  %  vac 

135  window.Et  =  cookiecut  (Et  ,D)  ;  %  turb 

136  %  collimate  field  at  tele  ap  (mult  by  conj  phase  from  mid  layer) 

137  Ev-col  =  window-Ev.  *  conj  ( Ev.mid.ph  ) ;  %  vac 

138  Et_col  =  window_Et .  *  conj  ( Ev_mid_ph  ) ;  %turb 

139  %  make  tele  ap  field  double  resolution  ,  convert  back  to  field 

140  Ev_col2  =  interp2(xt,  yt  ,  Ev_col.data  ,  xt2  ,  yt2  )  ; 

141  Ev_col2  =  emfield  (  Ev_col2  ,  wvl,  D)  ;  %  vac 

142  Et_col2  =  interp2(xt,  yt,  Et_col.data,  xt2,  yt2); 

143  Et_col2  =  emfield  (  Et_col2  ,  wvl,  D)  ;  %  turb 

144  %  propagate  fields  through  lenslets  to  detector  plane 

145  Ev_CCD  =  V-Llets  (  Llets  (  Ev_col2  )  ) ;  %  vac  field  at  detectors 

146  Et_CCD  =  V-Llets  (  Llets  (  Et_col2  )) ;  %  turb  field  at  detectors 

147  Iv.CCD  =  i  n  t  e  n  s  i  t  y  (Ev.CCD  )  ;  %  intensity  only  (vac) 

148  It_CCD  =  i  n  t  e  n  s  i  t  y  ( Et_CCD  )  ;  %  intensity  only  (turb) 

149 

150  9^%%%  manipulate  fields  in  Na  layer  (ground  to  sky  prop) 

151  It_Na  =  Iarray_Na_t(:,:,idx);  %  simple  rename  for  bookkeeping 

152  %  find  location  of  max  intensity  (centroid) 

153  yloc  =  -  round  (sum(sum(  It_Na).=!;([-N/2  :  N/2-1])... 

154  /sum(sum(  It_Na  )  ) ) )  ;  %  X  separation 

155  xloc  =  -  round  (sum(sum(  It_Na,2).*transpose([-N/2  :  N/2-1])... 

156  /sum(sum(  It_Na  )  ) ) )  ;  %  y  separation 

157  It_Na  =  circshift  ( It_Na  ,  [  xloc  yloc  ]  )  ;  %  center  max  intensity 

158  It_Na  =  emfield  ( It_Na  ,  wvl,  L);  %  convert  to  field  ,  irrad  already 

159  Iv.Na  =  emfield  ( Iarray_Na_v  ( :  , :  ,  idx  )  ,  wvl  ,  L);  %  irradiance  only 
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%cut  out  2mx2m  window  of  the  field 
window_Iv_Na  =  cookiecut  ( Iv_Na  ,D)  ;  %  vac 
window_It_Na  =  cookiecut  ( It _Na  ,D )  ;  %  turb 
%  demagnify  fields  from  Na  layer 

Iv_Na_dmag  =  magnify  ( window_Iv_Na  ,  0.025  );  %  vac 

It_Na_dmag  =  magnify  ( window_It_Na  ,  0.025  );  %  turb 

%  interpolate  demagnified  field  to  14  x  14  pixels 
Iv_Na_conv  =  interp2(xa,  ya ,  Iv_Na_dmag.data  ,  xb ,  yb )  ;  %  vac 

It_Na_conv  =  interp2(xa,  ya ,  It_Na_dmag.data  ,  xb ,  yb ) ;  %  turb 

98®o%%  convolution  to  make  elongated  spots  %%%%% 

temp_hor_v2  =  zer  os  ( 1  , 1 4 )  ' ;  %  initialize  horizontal  temp  variable 
temp_hor_t2  =  zeros  (1  ,14)  '; 

temp_ver_v2  =  zeros  ( 1  , 5 6  1 ) ;  %  initialize  vertical  temp  variable 
temp_ver_t2  =  zeros  ( 1  , 5 6  1 ) ; 

for  idx2  =  1:40;  %  horizontal  index  for  40  14x14  pixel  subaps 
for  idx3  =  1:40;  %  vertical  index  for  40  subaps 
%  pull  out  each  subap  from  array 
temp_hor_v  =  .  .  . 

Iv_CCD.data(idx*14-13:idx*14,idx2*14-13:idx2*14); 
temp_hor_t  =  .  .  . 

It_CCD.data(idx*14-13:idx*14,idx2*14-13:idx2*14); 

%  find  max  value  of  each  subap 
max_temp_hor_v  =  max(max(  temp_hor_v  ) )  ;  %  vac 
max_temp_hor_t  =  max(max(  temp_hor_t  ) ) ;  %  turb 
%  normalize  each  PSF  =  PDF  (divide  by  area  in  each  subap); 
temp_hor_v_norm  =  temp_hor_v/sum(sum(  temp_hor_v  )  )  ;  %  vac 
temp_hor_t_norm  =  temp_hor_t /sum(sum(  temp_hor_t  )  )  ;  %  turb 
%  convolve  PSF  with  pattern  in  sodium  layer 
temp_hor_v_c  =  max_temp_hor_v  .  .  . 

*convn  ( temp_hor_v_norm  ,  Iv_Na_conv  ,  '  same  '  ) ; 
temp_hor_t_c  =  max_temp_hor_t  .  .  . 

*  convn  (temp_hor_t_norm  ,It_Na_conv  ,  '  same  '  ) ; 

%store  convolved  patterns  horizontally 
temp_hor_v2  =  [  temp_hor_v2  ,  temp_hor_v_c  ]  ;  %  vac 
temp_hor_t2  =  [temp_hor_t2  ,  temp_hor_t_c  ]  ;  %  turb 

end 

%  store  each  row  of  subaps  vertically 
temp_ver_v2  =  [  temp_ver_v2  ;  temp_hor_v2  ]  ;  %  vac 
temp_ver_t2  =  [  temp_ver_t2  ;  temp_hor_t2  ]  ;  %  turb 
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temp_hor_v2  =  zeros  (1  ,14)  %  reset  variable  (vac) 

temp_hor_t2  =  z  er  os  ( 1  , 1 4 )  ' ;  %  reset  variable  (turb) 

end 

%  remove  initialized  zeros 
temp_ver_v2  =...  %  vac 

temp_ver_v2  (2:  length(temp_ver_v2)  ,2:  length(temp_ver_v2  )); 
temp_ver_t2  =...  %  turb 

temp_ver_t2  (2:  length(temp_ver_t2)  ,2:  length(temp_ver_t2  )); 
%  add  the  spot  locations 

Iv.final  =  Iv.final  +  temp_ver_v2  ;  %  vac 
It_final  =  It_final  +  temp_ver_t2  ;  %  turb 

end 

end 
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